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I We study, using transfer- matrix methods, the partition- function zeros of the 

square-lattice g-state Potts antiferromagnet at zero temperature (= square-lattice 
chromatic polynomial) for the boundary conditions that are obtained from an m x n 
grid with free boundary conditions by adjoining one new vertex adjacent to all the 
^ ' sites in the leftmost column and a second new vertex adjacent to all the sites in the 

. rightmost column. We provide numerical evidence that the partition-function zeros 

are becoming dense everywhere in the complex q-plane outside the limiting curve 
-SooCsq) for this model with ordinary (e.g. free or cylindrical) boundary conditions. 
Despite this, the infinite-volume free energy is perfectly analytic in this region. 
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1 Introduction 



It is well known that phase transitions do not occur in statistical-mechanical systems in 
finite volume, but occur only in the infinite- volume limit. One approach to studying phase 
transitions was introduced by Yang and Lee [72] in 1952, and involves investigating the 
zeros of the finite-volume partition function when one or more physical parameters (e.g. 
temperature or magnetic field) are allowed to take complex values. Lee and Yang showed, 
under mild conditions, that the accumulation points of these zeros in the infinite- volume 
limit constitute the only possible loci of phase transitions; away from such accumulation 
points of zeros, the infinite-volume free energy is analytic. 

We would like to stress the word "possible" in the preceding sentence. The finite- 
volume free energy is of course singular at any zero of the finite-volume partition func- 
tion (since logO — — oo); but it can happen that such a singularity disappears in the 
infinite- volume limit. To take a trivial example, suppose that the partition function 
in volume n is Zn{x) = x, where x is a parameter; then the finite- volume free energy 
fn{x) = log Zn{x) = n^^ logx is singular at a; = 0, but the infinite- volume free energy 
f{x) = lim fn{x) = is analytic at ,x = 0."*^ 

In this paper we would like to give an extreme example of this phenomenon, in which 
the zeros of the finite-volume partition functions appear to be dense in a large region 
of the complex plane, while the infinite-volume free energy is perfectly analytic there.^ 
Our example consists of the square-lattice Potts antiferromagnet at zero temperature — 
a model that we have studied in a series of previous papers [34-36,38,53,54] — but with 
some unusual boundary conditions. We begin by reviewing the needed background on 
Potts models and chromatic polynomials. 

The Potts model [47, 70, 71] on a regular lattice C is characterized by two parameters: 
the number q of Potts spin states, and the nearest-neighbor coupling v = e^"^ — l."^ Initially 
g is a positive integer and f is a real number in the interval [—1, +oo), but the Fortuin- 
Kasteleyn representation (reviewed in Section 2.1 below) shows that the partition function 
Zoiq, v) of the g'-state Potts model on any finite graph G is in fact a polynomial in q and 
V. This allows us to interpret q and v as taking arbitrary real or even complex values, 
and to study the phase diagram of the Potts model in the real (g, i')-plane or in complex 



One could object that in this case we actually have 




for X 7^ 
— oo for X = 



which is not analytic at x = (though the singularity at x = is removable). Here is a variant that 
avoids this objection: Z„{x) = x — 1/n. Then the finite- volume free energy /n(x) = n~^logZ„(x) is 
singular at x = 1/n, but it is easy to see that f{x) = lim fnix) = for all x (including x = 0). 

n— >oo 

^See also Section 8.2 for a brief discussion of the general question: When does an accumulation of 
partition-function zeros signal a nonanalyticity of the infinite-volume free energy? 

"^Here we are considering only the isotropic model, in which each nearest-neighbor edge is assigned the 
same coupling v. In a more refined analysis, one could put (for example) different couplings Vi,V2 on the 
horizontal and vertical edges of the square lattice, different couplings vi,V2,vs on the three orientations 
of edges of the triangular or hexagonal lattice, etc. 
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(g, f)-space. In particular, by studying the zeros of ZG{q,v) in complex (g,f)-space for 
larger and larger pieces of the lattice C, we can locate the possible loci of phase transitions 
in the real {q, f )-plane and more generally in complex (g, v)-spa.ce. 

The partition function for m x n lattices can be efficiently computed using transfer 
matrices. Though the dimension of the transfer matrix (and thus the computational 
complexity) grows exponentially in the width m — thereby restricting us in practice to 
widths m ^ 10-30 — it is straightforward, by iterating the transfer matrix, to handle quite 
large lengths n. Indeed, by implementing the transfer- matrix method symbolically (i.e., 
as polynomials in q and/or v) and using the Beraha-Kahane-Weiss theorem (reviewed in 
Section 2.3). wc can handle directly the limit n — > oo and compute the limiting curves 
Bm of partition-function zeros. At a second stage we attempt to extrapolate these curves 
to m = oo. 

Since the problem of computing the phase diagram in complex (g, w)-space is difficult, 
it has proven convenient to study first certain "slices" through (g, v)-space, in which one 
parameter is fixed (usually at a real vahic) while the remaining parameter is allowed to 
vary in the complex plane. One very interesting special case is the chromatic polynomial 
[v = —1), which corresponds to the zero-temperature limit of the Potts antiferromagnet 
{(3 J — — oo). In previous papers [34-36,38,53,54] we have used symbolic transfer-matrix 
methods to study the square-lattice and triangular-lattice chromatic polynomials for free, 
cylindrical, cyclic and toroidal boundary conditions. ^'^ Here we shall study the square- 
lattice model for the boundary conditions that are obtained from an rnp x uy square grid 
(m columns, n rows, free boundary conditions in both directions) by adjoining one new 
vertex adjacent to all the sites in the leftmost column and a second new vertex adjacent 
to all the sites in the rightmost column (see Figure la). 

The motivation for studying this family of graphs comes in part from work of one of 
the authors [63] on the generalized theta graphs [12]: namely, B*^**'^-* is defined to be the 
graph consisting of p chains in parallel between a pair of endvertices, each chain consisting 
of s edges in series (Figure 2). For this family, the following result holds: 

Theorem 1.1 [63, Theorems 1.1, 1.2 and 1.4] The set of chromatic roots of all the 
graphs Q^^'P^ is dense in the entire complex q-plane with the possible exception of the disc 
|g- 1| < 1. 

The proof of Theorem 1.1 proceeds in two steps: First, one fixes s, and shows that 
as p ^ oo the chromatic roots of the graphs Q^^'P^ accumulate densely on an explicitly 
computable real algebraic curve Cs- This argument is a straightforward application of the 
Beraha-Kahane-Weiss theorem. Second, one shows that as s — )■ oo the curves Cg become 
dense inC\{|g — 1|<1}. This requires an ad hoc (but not terribly difficult) argument. 
See Section 2.5 for a slightly more detailed summary. 

^See also the bibliographies of [34-36,38,53,54] for reference to the important related works of Shrock 
and collaborators. 

^We adopt Shrock's [56] terminology for boundary conditions: free (m-p x np), cylindrical (mp x np), 
cyclic (mp x np), toroidal (mp x np), Mobius (mp x nxp) and Klein bottle (mp x nxp). Here the 
first dimension (m) corresponds to the transverse ("short") direction, while the second dimension (n) 
corresponds to the longitudinal ( "long" ) direction. The subscripts F, P and TP denote free, periodic and 
twisted-periodic boundary conditions, respectively. 
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The computations arising in the proof of Theorem 1.1 are based on the series and 
parallel reduction laws for the Potts-model partition function [63,64]. This approach 
works because the graphs O^*'^') are series-parallel: we can compute the partition function 
ZQ(s,p){q,v) by computing the "effective coupling" v^s arising from putting s edges (each 
of weight v) in series and then placing p such chains in parallel. 

No such simple approach can work for the graphs Sm^n, because they are not series- 
parallel. Still, the graphs S^^n and the generalized theta graphs are closely related: 
indeed, the generalized theta graph 0("*+i'"') is obtained from Sm,n simply by deleting 
all the vertical edges. So it would not be surprising if the density-of-roots phenomenon 
observed for the generalized theta graphs occurred also for the family Sm,n- In this 
paper we would like to present some convincing evidence — which, however, falls short at 
present of a rigorous mathematical proof — that this is indeed the case (see in particular 
Conjecture 7.2 below). 

As in the analysis [63] of the generalized theta graphs, we proceed in two steps. First 
we fix m, and show that as n ^ oo the chromatic roots of the graphs Sm^n accumulate 
densely on a real algebraic curve Bm (as well as on certain isolated points). This argument 
is once again a straightforward application of the Beraha-Kahane- Weiss theorem. The 
more difficult problem is to compute the curves Bm and to study their behavior as m — )■ oo. 
In this paper we shall use the transfer- matrix formalism developed in [11,34,38,53,54] to 
compute the curves Bm for m < 6 (sec Figures 9-14 in Sections 6 and 7 below). We will 
see that the curve Bm contains 2m outward branches running to infinity in the complex 
g-plane and equally spaced in asymptotic angle. Moreover, these branches appear to 
become dense as m — > oo everywhere in the complex g-plane outside the limiting curve 
B(x,(sq) for this model with ordinary (e.g. free or cylindrical) boundary conditions — a 
curve that we have estimated numerically in previous papers [34,53] and which is entirely 
contained inside the disc Iq' — 1| < 3 (to give a crude bound). If this is indeed the case, 
it follows that the zeros of the chromatic polynomials of the graphs Sra,n are likewise 
becoming dense everywhere outside Soo(sq) and in particular everywhere in the region 
|g-l] > 3. 

As will be explained in Section 7 below (drawing on results from Section 4), these 
outward branches arise from the fact that the transfer matrix T"(m) for this model has 
two eigcnvahics — coming from two sectors in the block-diagonalization of T"(m) — that 
are almost equal at large q, with a relative difference of order q^"^, and which compete for 
dominance: the outward branches of the curves Bm, are nothing other than the loci where 
these two eigenvalues are equal in modulus. On the other hand, the difference between 
these two eigenvalues disappears as m — )■ cxd when |g| is large enough, with the result that 
the infinite- volume free energy is perfectly analytic at large 

A behavior of this type was observed recently by Jacobsen, Richard and Salas [33] for 
the Potts model with cyclic boundary conditions (i.e. periodic boundary conditions in the 
longitudinal direction) when q is fixed at a Beraha number and the temperature variable 
V is taken complex. On the square (resp. triangular) lattice, there are two eigenvalues — 
coming once again from two sectors in the block-diagonalization of the transfer matrix — 
that are almost equal at large v, with a relative difference of order f ^™ (resp. v"^'^"^~^^). 
This leads to 2m (resp. 4m — 2) outward branches in the curve B^. that appear to become 
dense as m — > oo everywhere outside a bounded region in the complex v-plane. These 
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same outward branches were observed earlier by Shrock and Chang for m = 1,2 [55] and 
m = 3 [15] on the square lattice and for m = 2 on the triangular lattice [14]; they exist 
for all q ^ 0,1, not only Beraha numbers. 

The principle underlying all these examples is a general one: Suppose that the transfer 
matrix at some parameter value xq has two dominant eigenvalues of equal moduhis, call 
them Ai(xo) and \2{xq) = e*^Ai(a;o), and that for parameters x ^ xq we have X2{x) — 
e^^Xi{x) = A{x — Xq)'^ + o{{x — Xq)'^) with ^4 7^ and some integer k > 1 [with the obvious 
modifications in case Xq = 00]. Then the curve contains 2k branches passing through 
Xq, with an asymptotic angle 7r/k between neighboring branches. If A; = 1 this is simply 
the generic situation for equimodular curves; if A; > 1 it corresponds to a higher-order 
crossing. In particular, if k tends to infinity as m — )■ 00, these branches are likely to 
become dense in a neighborhood of xq. 

The plan of this paper is as follows: In Section 2 we review some needed background 
on chromatic and Tutte polynomials, transfer matrices, and the Beraha-Kahane-Weiss 
theorem. In Section 3 we present the transfer-matrix theory for the family Sm,n and we 
prove some properties of the dominant diagonal elements and the dominant eigenvalues. 
In Section 4 we present the results of our transfer-matrix computation of the large-g 
expansion of the leading eigenvalues for 1 < m < 11 and note some surprising properties. 
In Section 5 we study the thermodynamic limit m — > 00 of the corresponding free energies. 
In Section 6 wc present the numerically computed limiting curves Bm for widths 1 < m < 
6. In Section 7 we present some conjectures on the behavior of these limiting curves as 
m — )■ 00 and we relate these conjectures to the properties of the eigenvalues that were 
empirically observed in Section 4. Finally, in Section 8 we conclude with some brief 
discussion and suggestions for future research. In Appendix A we prove Proposition 3.5 
concerning the dominant diagonal entries in the transfer matrix. In Appendix B we 
compute a general formula for the dimensions of our transfer matrices as a function of 
the width m. In Appendix C we present the transfer-matrix theory for the family Sm,n 
obtained by a 90° rotation of the grid; this plays an important role in checking the 
correctness of our transfer-matrix computations. 

2 Preliminaries 

In this section we review briefly some needed background on chromatic and Tutte 
polynomials (Section 2.1), transfer matrices (Section 2.2), and the Beraha-Kahane-Weiss 
theorem (Section 2.3). We also recall the easy analysis of the chromatic roots of bi-fans 
and bipyramids (Section 2.4) and briefly summarize the argument from [63] concerning 
the chromatic roots of generahzed theta graphs (Section 2.5); both of these will serve as 
models for our analysis here of the graphs S^^n- 
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2.1 Chromatic polynomials, Potts models, and all that 



Let G = {V, E) be a finite undirected graph, and let g be a positive integer. Then the 
q-state Potts-model partition function for the graph G is defined by the Hamiltonian 

^^Potts(o-) = - Je5{ai,aj) , (2.1) 

e=ij&E 

where the spins cr = {ai}i^v take values in {1,2, ... ,q}, the Jg are coupling constants, 
and the 5 is the Kronecker delta 

r/ ,N I 1 if a = 6 _ ^, 

The partition function can then be written as 

^G°"^(?,v) = Yl n + ' (2-3) 

a: V^{l,2,...,q} e=ijeE 

where Ve = e^'^" — 1. Please note, in particular, that if we set Ve = —1 for all edges e, 
then Zq"^^'^ gives weight 1 to each proper coloring and weight to each improper coloring, 
and so counts the proper colorings. Proper g-colorings {Vg = —1) thus correspond to the 
zero-temperature {(3 +00) limit of the antiferromagnetic ( Jg < 0) Potts model. 

It is far from obvious that ZQ°*^^{q,v), which is defined separately for each positive 
integer q, is in fact the restriction to e Z+ of a polynomial in q. But this is in fact the 
case, and indeed we have: 

Theorem 2.1 (Fortuin— Kasteleyn [27,42] representation of the Potts model) 

For integer q > 1, 

^G°"1?,v) = 5^ g^^(-^) n (2.4) 

ACE eeA 

where k{A) denotes the number of connected components in the subgraph {V,A). 



Proof. In (2.3), expand out the product over e & E, and let A C E be the set of edges for 
which the term VeS{ai, aj) is taken. Now perform the sum over maps a: V ^ {1,2, . . . ,q}: 
in each component of the subgraph [V, A) the color cTj must be constant, and there are 
no other constraints. We immediately obtain (2.4). □ 

Historical Remark. The subgraph expansion (2.4) was discovered by Birkhoff [10] 
and Whitney [69] for the special case Ve = —1 (see also Tutte [67,68]); in its general form 
it is due to Fortuin and Kasteleyn [27,42] (see also [21]). 

The foregoing considerations motivate defining the multivariate Tutte polynomial of 
the graph G: 

^G(g,v) = E^'^^^n^- (2.5) 

ACE eeA 
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where q and v = {ve}eeE are commuting indeterminatcs. If we set all the edge weights 
Ve equal to the same value we obtain a two-variable polynomial that is equivalent to 
the standard Tutte polynomial Tq^x^h) after a simple change of variables (see [64]). If 
we set all the edge weights equal to —1, we obtain the chromatic polynomial Pg{(i) — 



Further information on the multivariate Tutte polynomial Zoiq, v) can be found in a 
recent survey article [64]. 

2.2 Transfer matrices 

For any family of graphs G„ = (Ki, En) consisting of n identical "layers" with identical 
connections between adjacent layers, the multivariate Tutte polynomials of the Gn (with 
edge weights likewise repeated from layer to layer) can be written in terms of a transfer 
matrix [11,53]. Here we briefly summarize the needed formalism [53] specialized to the 
case of an m X n square lattice with free boundary conditions in both directions. The 
modifications needed to handle the extra sites at left and right in the graphs Sm^n will be 
discussed in Section 3. 

Consider the m x n square grid with edge weights Vi^i+i on the horizontal edges (1 < 
i < m — 1) and Vi on the vertical edges {1 < i < m). We fix the "width" m and consider 
the family of graphs Gn obtained by varying the "length" n; our goal is to calculate the 
multivariate Tutte polynomials ZG!^(g, v) for this family by building up the graph Gn 
layer by layer. What makes this a bit tricky is the nonlocality of the factor g^'(^) in (2.5). 
At the end we will need to know the number of connected components in the subgraph 
(Vn, A); in order to be able to compute this, we shall keep track, as we go along, of which 
sites in the current "top" layer are connected to which other sites in that layer by a path 
of occupied edges (i.e. edges of A) in lower layers. Thus, we shall work in the basis of 
connectivities of the top layer, whose basis elements e-p are indexed by partitions V of 
the single-layer vertex set {1, . . . , m}. The elementary operators we shall need are: 

• The join operators 



where V»ij is the partition obtained from V by amalgamating the blocks containing 
i and j (if they were not already in the same block) . Note that all these operators 
commute. 

• The detach operators 



where V\i is the partition obtained from V by detaching i from its block (and thus 
making it a singleton) . Note that these operators commute as well. 

Note, finally, that commutes with Jjj whenever k ^ {i,j}- 

The horizontal transfer matrix, which adds a row of horizontal edges, is 

m— 1 



(2.6) 




(2.7) 



H = JJ(1 + Vj,j+iJj,j+i) 



(2.8) 



i=l 
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(note that all the operators in the product commute). The vertical transfer matrix, which 
adds a new row of sites along with the corresponding vertical edges, is 



y^HivJ + Di) (2.9) 



1=1 



(note once again that all the operators commute). The multivariate Tutte polynomial for 
Gn is then given [53] by the formula 

ZG^{q,v) = a;^H(VH)"-ieid , (2.10) 

where "id" denotes the partition in which each site ie{l,...,m}isa singleton, and the 
"end vector" u?'^ is defined by 

= g"'' ■ (2.11) 

The transfer matrix is thus 

T = VH . (2.12) 

In what follows we shall use a convenient shorthand notation for the basis vectors e-p: 
namely, we denote by 1 the basis element Cid corresponding to the partition in which each 
site is a singleton, and we denote the action on eid of a join operator Jjj by the Kronecker 
deltas Sij. Thus, for instance, for m = 3 we have the five basis vectors 

eid = e{{i},{2},{3}} = 1 (2.13a) 

e{{i,2},{3}} = Si2 (2.13b) 

e{{i,3},{2}} = Si3 (2.13c) 

e{{i},{2,3}} = S23 (2.13d) 

e{ {1,2,3}} = ^123 = 5l25l3523 = C)l2^13 = f'l2(^"23 = f'l3<523 (2.13e) 

For m > 4 we also have basis vectors that are products of deltas, such as 513^24 and so 
forth. This notation allows describing partitions and their corresponding basis vectors in 
a fairly compact way. 

In principle we are working here in the space spanned by the basis vectors e-p for 
all partitions V of {1, . . . ,m}; the dimension of this space is given by the Bell number 
Bm [18,61,65,66]. However, it is easy to see, on topological grounds (thanks to the 
planarity of the Gn), that only non- crossing partitions can arise. (A partition is said to 
be non-crossing ifa<b<c<d with a, c in the same block and b, d in the same block 
imply that a, b, c, d are all in the same block.) The number of non-crossing partitions of 
{1, . . . , m} is given by the Catalan number [61, 66] 

C„ . . ^ f^'") . (2.14) 



m!(m + l)! m+l\m 

When the horizontal couplings are all equal to —1 (which is the case for the 

chromatic polynomial), then the horizontal operator H is a projection, and we can work 
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in its image subspace by using the modified transfer matrix T' — HVH in place of T = VH, 
and using tlie basis vectors 

fp = Hep (2.15) 

in place of e-p. Please note that fp = if P has any pair of nearest neighbors in the 
same block. We thus work in the space spanned by the basis vectors fp where is a 
non-crossing non- nearest-neighbor partition of {1, . . . , m}. The dimension of this space is 
given [60, Proposition 3.6] [45] by the Motzkin number M^_i, where [9,19,29,51,61,66]^ 



Ln/2J 

= E 

3=0 



(2.16) 



Finally, spatial symmetries further restrict the subspace whenever the couplings ^1,1+1 
and Vi are invariant under the symmetry. Here the relevant symmetry is reflection in 
the center line of the strip. For reflection-invariant couplings, we can work in the space 
of equivalence classes of non-crossing non-ncarcst-ncighbor partitions modulo reflection. 
We denote the corresponding symmetry-reduced transfer matrix by T", and its dimension 
SqFree(m) is then given by the number of these equivalence classes. The exact expression 
for SqFree(m) was obtained in [13, Theorem 2]: 



SqFree(m) = ]-M^_i + 



where 



(m'- 1): 



m + 1 



\m'/2\ 
j=0 



m — J 



(j!)2(m'-2j)! 



m 



(2.17) 



(2.18) 



and [pj stands for the largest integer < p. (We give a new proof of this formula in 
Appendix B.) The asymptotic behavior of SqFree(m) is given by [13, Corollary 1] 



SqFree(m) 



3'^m 



-3/2 



1 + 



as m — >■ 00 



(2.19) 



The values of all these dimensions for m < 16 are displayed in Table 1. 

Some further structural properties of the transfer matrices can be found in Ref . [54] . 



2.3 Beraha—Kahane— Weiss theorem 

A central role in our work is played by a theorem on analytic functions due to Bc- 
raha, Kahane and Weiss [5-8] and generalized slightly by one of us [63]. The situation 
is as follows: Let D be a domain (connected open set) in the complex plane, and let 
«!, . . . , QiM, Ai, . . . , Am {M > 2) be analytic functions on D, none of which is identically 
zero. For each integer n > 0, define 

M 

fn{z) = ^afe(^)Afe(^)V (2.20) 
fc=i 

^ Warning: Several references use the notation nin to denote what we call M„; and one reference [19] 
writes M„ to denote a different sequence. 
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We are interested in the zero sets 

Z{U) ^ {zeD: fn{z)^0} (2.21) 

and in particular in their hmit sets as n ^ oo: 

hminf = {z & D: every neighborhood U 3 z has a nonempty intersection 

with all but finitely many of the sets (2.22) 

limsup.E(/„) = {z E D: every neighborhood U 3 z has a nonempty intersection 

with infinitely many of the sets Z(fn)} (2.23) 

Let us call an index k dominant at z if |Ajt(2;)| > |A;(2;)| for alH (1 < i < M); and let us 
write 

-Dfc = {z E D: k is dominant at z} . (2.24) 
Then the limiting zero sets can be completely characterized as follows: 

Theorem 2.2 [5—8,63] Let D be a domain in C, and let ai, . . . , au, Ai, . . . , Am (M > 

2) be analytic functions on D, none of which is identically zero. Let us further assume 
a "no- degenerate- dominance" condition: there do not exist indices k ^ k' such that \k = 
cuXk' for some constant cu with = 1 and such that (= D^') has nonempty interior. 
For each integer n > 0, define fn by 

M 

fniz) = 5]afe(z)Afe(z)V 

k=l 

Then liminf Z{fn) — limsup 2^(/„), and a point z lies in this set if and only if either 

(a) There is a unique dominant index k at z, and ak{z) = 0; or 

(b) There are two or more dominant indices at z. 

Note that case (a) consists of isolated points in D; we refer to these as isolated limiting 
points. Case (b) consists of curves (plus possibly isolated points where all the A^ vanish si- 
multaneously); we call these limiting curves or dominant equimodular curves. Henceforth 
we shall denote by B the locus of points satisfying condition (b). 

We shall often refer to the functions A^ as "eigenvalues", because that is how they 
arise in the transfer-matrix formalism. 

2.4 Chromatic roots of bi-fans and bipyramids 

The simplest case of the graphs Sm,n is width m = 1: the graphs 5*1^^ are called bi-fans 
(Figure 3a), and their chromatic polynomials can be easily computed. Here we would like 
to recall this computation and to discuss the behavior of the chromatic roots as n — >■ oo. 
Similar methods apply to a closely related class of graphs called bipyramids (Figure 3b). 

We recall that the join G + H oi two graphs is the graph obtained from the disjoint 
union of G and H by adding one edge connecting each pair of vertices i e ViG)., j G V{H). 
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Let K2 be the graph with two vertices and no edges, and let G be any graph. Then 
the chromatic polynomial of the join G + K2 is given in terms of the chromatic polynomial 
of G by the easy formula 

Pg+kM) = (iPaiq - 1) + q{q - ^)PG{q - 2) . (2.25) 

Indeed, the first (resp. second) term on the right-hand side of (2.25) counts the proper 
g-colorings of G + K2 in which the two vertices of K2 are colored the same (resp. different). 

Using the well-known chromatic polynomials of the n-vertex path P„ and the n-vertex 
cycle Cn, 

PpM) = (2.26) 

PcM) = (9-ir + (9-i)(-ir (2.27) 

we obtain the chromatic polynomials of the bi-fan P„ -|- K2 and the bipyramid Gn + K2: 

Pp.+kM = g(g-l)(g-2)"-^ + g(g-l)(g-2)(g-3)"-^ (2.28) 

Pc„+kM = ?(?-2)" + g(g-l)(g-3)" + g(g'-3g + l)(-l)'^ (2.29) 

Using the Beraha-Kahane-Weiss theorem, we see immediately that the chromatic 
roots of the bi-fans accumulate as — )■ cxd on the curve \q — 2\ = \q — 3\, i.e. the vertical line 
Reg = 5/2 [58]. (There are also isolated limiting points at g = 0, 1,2.) The analysis for 
the bipyramids is slightly more complicated, because there are three "eigenvalues" {q — 2, 
g — 3 and —1) that need to be compared. Elementary computations show [50, 57] that 
the limiting curve B consists of the pair of semi-infinite vertical line segments Reg = 5/2, 
I Img| > v^/2 together with the pair of circular arcs g = 2 -|- e'^ and g = 3 — e*^ for 
-7r/3 <e< 7r/3.^ 

It is worth remarking that if one plots the chromatic roots of the bi-fans or bipyramids, 
one finds a large number of roots accumulating on the parts of B fairly close to the real 
axis, together with a few roots at large imaginary part lying far to the right of the line 
Reg = 5/2.^ (In Figure 4 we depict the chromatic roots of the n = 100 bi-fan. The 
rightmost zeros are located at g ^ 21.8463758923 ± 14.9443992222 «.) This behavior 
is at first surprising (at least it was to us), but after some thought one realizes that 
it is perfectly consistent with the Beraha-Kahane-Weiss theorem: the hmiting curve 
Reg = 5/2 contains the point at infinity, and roots can tend to infinity in the topology of 
the Riemann sphere (which turns out to be the relevant sense) in many different ways; in 
particular, their real parts need not tend to 5/2. In fact, it turns out that the rightmost 
root has 

n 



logn — log log n 



\ ^ 7^^ ^ ^/^loglogn 



logn \ log n 



(2.30) 



''Read and Royle's formula for the chromatic polynomial of the bipyramids [50, p. 1017] has a typo- 
graphical error: — 3z + 3 should be z'^ — 3z + 1. 

^For the bipyramids this was noticed already two decades ago by Read and Royle [50, pp. 1019-1020]. 
It was further observed by Shrock and Tsai [58, Sections III. A and III.B] for some classes of graphs that 
generalize the bi-fans and bipyramids. 
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so that its real and imaginary parts both tend to infinity as n — )■ oo. This asymptotic 
formula can be proven by a minor modification of the analysis given in [12, Section 5] for 
the chromatic roots of the complete bipartite graphs i^2,n — -^n + -^2- 

2.5 Chromatic roots of generalized theta graphs 

Let us now summarize briefly the argument from [63] concerning the chromatic roots 
of the generalized theta graphs. 

The chromatic polynomial of is [12,59,63] 

Pei^AQ) = Q>^=AQy + QiQ ' (2.31) 

where 

A,,.(,) = (^-^r-(-^r (2.33) 

Indeed, \=,s{q) [resp. A^^s(g)] is the number of ways of coloring the internal vertices of 
one s-edge path when the two endvertices are colored the same (resp. different). Please 
note that the leading large-g term in both X=,s{q) and Xjt,s{q) is q^"^ = q^, and that they 
start to differ at order q*^: indeed, we have the curious formula 

X=,s{q) - X^Aq) = (-1)^' • (2-34) 

Let us begin by fixing s and considering the limit p — )■ oo. Using the Beraha-Kahane- 
Weiss theorem, we see immediately that the chromatic roots of the graphs Q^^'P^ accu- 
mulate as p — )■ oo on the real algebraic curve curve Cg defined by |A=^5(g)| = |A^^5(g)|, or 
equivalently 

Cs-. \1 + {q-l){l-q)-'\ = |l-(l-g)-^| . (2.35) 
To handle the limit s — > oo, we use the following lemma: 

Lemma 2.3 [63, Lemma 1.6] Let Fi,F2,G be analytic functions on a disc \z\ < R 
satisfying \G{0)\ < 1 and G ^ constant. Then, for each e > 0, there exists sq < oo such 
that for all integers s > Sq the equation 

\l + F,{z)G{zy\ = |1 + F2(^)G(^)^| (2.36) 

has a solution in the disc \z\ < e. 

This lemma immediately implies that the curves Cs accumulate as s — )■ oo at every point in 
the region \q — l\ > 1. A strengthened version of this lemma is proven in [63, Theorem 4.1]. 

In Figure 5a we show the curve Cg together with the chromatic roots of for 
s = 5 and p = 25, 100,400. Just as for the bi-fans, we see clearly the accumulation of 
roots on the inner parts of the curve Cs, together with a few points at larger |g| lying far 
off this curve. All this can presumably be explained by a generalization of the analysis 
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given in [12, Section 5] for the special case s = 2. In particular, we expect that the root 
of largest modulus grows like {p/ logp)^^^^~^^ as p ^ oo.^ This explains the extremely 
slow convergence of the roots to the outer parts of the curve Cg in Figure 5a. Note 
in particular that a naive computational approach, without the insight provided by the 
Beraha-Kahane-Weiss theorem, would give httle indication that the roots are in fact 
unbounded as p — )■ oo. 

In Figure 5b we show the curves for s = 5, 10, 15. We see clearly how these curves 
are becoming dense in the entire complex plane with the exception of the disc |? — 1| < 1. 
The curves Cg do enter this disc, but as s — > oo they retreat to its boundary, so that no 
point of the open disc is a limit point of the curves Cg. 

Finally, let us note, for future reference, the formula for the multivariate Tutte poly- 
nomial of 0(^'P) when the same weights vi, . . . ,Vs (but not necessarily vi — . . . — Vg = —1) 
are assigned to each path [63]: 



where 



A=,s(g,v) 



We have the curious formula 



i=l 



s s 

n(^+^) -n 



.1=1 



i=l 



(2.37) 

(2.38) 
(2.39) 



A=,3(g,v) - A^,,(g,v) = n^'' 



which generalizes (2.34). 



i=l 



(2.40) 



3 Transfer-matrix theory for the family Sm,n 

In Section 2.2 we explained the transfer-matrix formalism for computing the multi- 
variate Tutte polynomial of an m x n square lattice with free boundary conditions in both 
directions. In this section we begin (Section 3.1) by discussing the modifications needed 
to handle the extra sites at left and right in the graphs Sm^n- We then (Section 3.2) go on 
to prove some structural properties of the transfer matrix for the family Sm,n, by analogy 
with what was done in Ref. [54] for free and cylindrical boundary conditions; in particular 
wc show that the transfer matrix T can be decomposed into two independent blocks T= 
and T^. Next (Section 3.3) we introduce the families of graphs S'^„ and which 
shed further light on the block-diagonalization of the transfer matrix for the family Sm,n- 
Then (Section 3.4) we prove that the large- |g| expansion for the dominant eigenvalue of 

^Brown, Hickman, Sokal and Wagner [12, Theorem 1.3] proved that the roots are bounded in modulus 
by const x p/ logp, uniformly in s and p. But this bound is presumably far from sharp when s > 2. 
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the transfer matrix T= or holds uniformly in the strip width m (throughout the an- 
tiferromagnetic regime), by invoking results from [62]. Finally (Section 3.5), we prove a 
key identity concerning the difference between the dominant eigenvalues A* -(m, q, v) and 
A*,^(m,g,v) corresponding to S=^^ and S^^^. 

3.1 Basic method 

We want to compute the chromatic polynomial Psm,niQ) — more generally, the 
multivariate Tutte polynomial Zs^„{q,v) with periodically repeated edge weights — for 
the graph Sm,n that is obtained from a square-lattice strip of width m and length n (with 
free boundary conditions in both directions) by connecting every site in the leftmost (resp. 
rightmost) column to an extra site s (resp. s'). Thus, the number of sites of this strip is 
mn + 2 (see Figure la). 

We shall also denote the strip Sm,n as mF++ x rip: here the subscript F denotes free 
boundary conditions, while the subscript F++ denotes free boundary conditions with 
the two extra sites attached at left and right. A square-lattice strip with standard free 
boundary conditions in both directions will be denoted as m-p x rip. 

In the strip Sjn,n we assign weight Vo^i to all the edges connecting the leftmost site 
(which we can think of as "column 0") to the sites of column 1, weights to the 

horizontal edges linking column i to column i + 1 {1 < i < m — 1), weight Vm,m+i to all 
the edges connecting the sites of column m to the rightmost site (which we can think of 
as "column m + 1"), and weight Vi to the vertical edges in column i {1 < i < m). Of 
course, for the chromatic polynomial all these weights will be taken to be —1. 

The computation of the multivariate Tutte polynomial for this strip family can be 
achieved by using a minor modification of the transfer-matrix method explained in Sec- 
tion 2.2. As before, we fix the "width" m and consider the family of graphs Sm,n obtained 
by varying the "length" n; our goal is to calculate the multivariate Tutte polynomials 
Zsm,„{Q,^) for this family by building up the graphs Sm,n layer by layer. 

One way of constructing the transfer matrix T(m) for the family 5"^ „ is to treat this 
family as a special case of a square-lattice strip {m + 2)p x rip: see Figure lb. Thus, we 
choose to work in the space of connectivities {e-p} associated to the non-crossing partitions 
V of the set {0,1,2, ... ,m,m + 1}. The dimension of this space is the Catalan number 



On the edges of the {m + 2)p x np strip that correspond to the edges of the starting 
graph Sm,n, we simply use the desired weights {0 < i < m) and Vi {1 < i < m). On 
the vertical edges in the first and last columns (depicted with dashed lines in Figure lb), 
we associate a weight Ve +oo, which corresponds to contracting all the sites in each 
of these columns into a single site. In other words, we use two different vertical-bond 
operators depending on which column we are in: 




(3.1) 



m 



H = JJ(1 Vj,j+iJj,j+i) 



(3.2) 



i=0 
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[compare (2.8)] and the vertical transfer matrix is given by 

m 

V = \[{v,I+Di). (3.3) 

j=i 

Note that the formula for V is identical to the formula (2.9) for an ordinary mp x np 
grid; the only difference is that here the operators act on the space of partitions of 
{0, 1, 2, . . . , m, m + 1} rather than {1,2...., m}. 

The full transfer matrix is then T = VH, and the multivariate Tutte polynomial 
Zsm,M^^) is given by 

Zs^,M.^) = a;THT'^-ieid , (3.4) 
where denotes the basis vector associated to the partition 

id = {{0},{l},...,{m},{m + l}} (3.5) 

and the left vector w acts on a connectivity state as in (2.11). 

When the horizontal couplings are all equal to —1 (as is the case for the chromatic 
polynomial), the horizontal operator H is a projection, so that (as explained in Section 2.2) 
we can work in its image subspace by using the modified transfer matrix T' = HVH and 
the basis vectors fp = Hep, where V runs over the non-crossing non-nearest-neighbor 
partitions of{0,...,m-|-l}. The dimension of this space is the Motzkin number Mj^+i 
[cf. (2.16)]. We then write the chromatic polynomial for Sjn,n as 

PsU^) = u:^{rr-'fi, . (3.6) 

Finally, when the couplings are also reflection-symmetric (e.g., when they are all — 1), 
we can work in the space of equivalence classes of non-crossing non-nearest-neighbor 
partitions modulo reflection. The dimension of this space is given by SqPree(m -|- 2) [cf. 
(2.17)]. Hereafter, we will denote by T"(m) this symmetry-reduced transfer matrix. 

3.2 Block-diagonalization of the transfer matrix (I) 

One important property of the transfer matrix T (or T' or T") for the family Sm,n is 
that it can be decomposed into two independent transfer matrices of a lower dimension. 
To sec this, it suffices to note that the projection Jo,m+i commutes with both H and V: it 
commutes with H because all join operators commute with each other, and it commutes 
with V because (3.3) does not contain the detach operators Dq or D^+i. Therefore, T 
(or T' or T") block-diagonahzes if we choose a basis in which the first group of basis 
elements corresponds to the range of Jo,m+i, and the second group corresponds to the 
range of / — Jo,m+i- We call these two blocks T= and T^, respectively. Indeed, this 
block-diagonalization holds for any transfer matrix built out of operators of the general 
form 

H = ^ CA]^Ji,i+i (3.7a) 

AC{0,...,m} isA 

V = E ^^n^^ (3.7b) 

BC{l,...,m} ieB 
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irrespective of the coefficients {ca} and {(is}- The point, once again, is simply that V 
does not contain the detach operators Dq or D^+i. 

Let us now recall that id denotes the partition of {0, . . . , m + 1} in which each element 
is a singleton, and that id • (0, m + 1) denotes the partition in which {0, m + 1} is a block 
and all the other elements are singletons. Let us call a partition V of {0, . . . ,m + 1} 
non-trivial if it is neither id nor id • (0, m + 1). 

Proposition 3.1 For any operators H and\/ of the form (3.7), we can write 

VHJo,rn+ieid = ^=Jo,m+ieid + 0'V^Q,m+l^V (3.8a) 

V non-trivial 

VH (7 - Jo,r„+i)eid = (-/" - Jo,m-n)eid + ^ (7 - Jo,mH-i)ep (3.8b) 

V non-trivial 

for some coefficients t=, t^, {a^} and {a^} that are polynomials in q, {c^} and {d^}. 
Furthermore, these coefficients satisfy the identity 

t= — t^ = C{o,...,m}C?{l,...,m} (3.9) 

reminiscent of (2.34)/ (240). 



Proof. Let us begin by considering the quantity Heid: 

Heid = X] Ji,i-Hi j eid . (3.10) 

Let us now apply the operator V to (3.10). It is clear that wc obtain three types of terms: 
a) terms that contain non-trivial partitions b) terms that contain the trivial partition 
id, and c) a single term that contains the trivial partition id • (0, m + 1). Indeed, this 
latter term can arise in only one way: first we must join all the sites {0, . . . , m + 1} into a 
single block (i.e., take A — {{),... , m} in H), and then we must detach all the inner sites 
(i.e., take B — {1, . . . , m} in V). It follows that 

VHCid = C{0,...,m}d{l,...,m}^0,m+l^id + OidGid + ^ a-p G-p (3.11) 

non-trivial 

for some coefficients Oid and {ap} that will clearly be polynomials in g, {ca} and {ds}- 
From this equation is easy to derive (3.8), as the join operator Jo,m-i-i commutes with all 
the join and detach operators arising in (3.11). Indeed, using Jo,r„+i = Jo,m-n we find 

t= = aid + C{o,...,m}C^{i,...,m} (3.12a) 
t^ = Oid (3.12b) 

Prom this relation, (3.9) follows immediately. □ 
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When H and V are the transfer matrices (3.2)/(3.3) for the multivariate Tutte poly- 
nomial, we have C|o,...,m} = Hl^Io'^i.^+i ^{^,--;m} = 1) so that (3.9) becomes 

m 

t=-u = • (3-13) 

Note that this formula holds irrespective of the values of the vertical couplings Vi. In par- 
ticular, it holds when Vi = 0; this explains the "curious formula" (2.40) for the generalized 
theta graphs. ^° When all the horizontal couplings I'j^j+i equal —1, we have 

t=-t^ = {-ir+\ (3.14) 

which includes (2.34) as a special case. 

The next result is similar to Proposition 3.1 in [54]. 

Proposition 3.2 Consider operators H and V of the form (3.7) where the coefficients 
{ca} and {ds} are numbers (i.e., independent of q). Then, using the notation of Propo- 
sition 3.1, 

t= = C0d^i^,,,^jnjq"^ + terms of order at most q"^^^ (3.15a) 
t^ — C0d{i^,„^rn}(i^ + terms of order at most q^~^ (3.15b) 

Furthermore, all the coefficients {a^} and {a^} defined in Proposition 3.1 are polynomials 
in q of degree at most m — 1. 



Proof. First of all, it is obvious that each entry in the transfer matrix T = VH is a 
polynomial in q. Indeed, from (2.7) it is clear that we get a factor of q every time we 
apply the operator Dj to a partition in which i is a singleton. Indeed, wc can maximize 
the number of factors of q by applying the vertical transfer matrix V to a partition in 



^^Thcrc is a small missing step here, because (3.13) refers to the dominant diagonal entries of the 
transfer matrix, while (2.40) refers to the transfer-matrix eigenvalues. But in the case of the generalized 
theta graphs Q'^^'P^ these two quantities coincide, as can easily be seen by expressing Zq(s,p) (g, v) using 
a transfer-matrix formalism. Labelling the vertices along each of the p chains from to s, it is not hard 
to see that the partition function can be written in the basis {1, 5o,s} as follows: 



.^0(.,p) (g, v) 



„2 \T / A^,s(g,v) \ / 1 
y \ I s-i 1 I ^ 



n A=,s(g,v) ] \^ 

i=0 



Since the transfer matrix is lower-triangular, the leading diagonal entries of the transfer matrix coincide 
with its eigenvalues. 

Let us remark that if we choose the basis {5os, 1 — <5os} = {Jos, {I — Jos)}eid, then the above expression 
simplifies to 

^e,..(„v,.,(,!OM"'t"' a,v))'(i 

and we obtain (2.37). 
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which every site 1, . . . ,m is a singleton, i.e. either eid or Jo,m+ieid- In particular, from 
(3.7b) we have 

BC{l,...,m} 

— {d{i,...,m}Q."^ + terms of order at most q"^~^) Bid (3.16b) 

and analogously for Jo,m-t-ieid- If we apply the vertical transfer matrix to any non-trivial 
partition, we get a polynomial in q of degree at most m — l. Now let M denote any linear 
combination of Jo,m+i a-nd / — Jo,m+i) a-nd consider 

HMeid = Yl (II-'m+i) (^-l^^) 

AC{0,...,m} \i&A / 

= Ce,Meid + Y bpMe-p (3.17b) 

■p non-trivial 

for some quantities b-p that are polynomials in {ca} (and of course independent of q). 
Using (3.16)7(3.17) it is obvious that 

VHMeid = C0 Y dBq\^^MeiA + J2^'p{q)e-p (3.18) 

BC{l,...,m} V 

where the coefficients b'pi^q) are polynomials in q of degree at most m — l. □ 

In view of Proposition 3.2, we shall henceforth refer to t= and as the "dominant 
diagonal terms" in the transfer matrix, as they are indeed dominant at large Further- 
more, we can deduce from Propositions 3.1 and 3.2 the leading large-|g| behavior of the 
eigenvalues. We begin with a simple perturbation lemma proved in Ref. [54]: 

Lemma 3.3 [54, Lemma 3.2] Consider an N x N matrix M{^) = {Mij{^))f^j^^ whose 
entries are analytic functions of ^ in some disc \^\ < R. Suppose that Mu = 1 and that 
Mij — 0(C) for {i,j) 7^ (1)1)- Then, in some disc \^\ < R' , M(C) has a simple leading 
eigenvalue //*(0 that is given by a convergent expansion 

oo 

/..(O = 1 + E«^^' (3-19) 

k=2 

[note that ai — 0] with associated eigenvector 

oo 

v.(0 = ei + J]v,e', (3.20) 

k=l 

while all other eigenvalues are 0{^). 

Remarks. 1. The key fact here is that the eigenvalue shift in (3.19) begins at order 
not order ^. 
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2. The "small" eigenvalues need not be analytic in ^. For instance, 




M(0 = H (3.21) 



has eigenvalues /i — l and /i — ±^^/^. □ 

To apply Lemma 3.3 to our transfer matrix T= (resp. T^), we set ^ = and 
M = T=/t= (resp. T^/t^). We then have:^^ 

Corollary 3.4 Consider operators H and V of the form (3.7) where the coefficients {ca} 
and {ds} are numbers (i.e., independent of q) and C0d^i^,,,^m} 0- ThenT= andJ^ each 
have a single eigenvalue that is analytic for large \q\ and behaves there like const x q"^: 
more precisely, these eigenvalues have convergent expansions 



A. 



t= 



t 



= 1 + (3.22a) 

k=2 
oo 

= 1 + (3-22b) 



k=2 



SO that in particular 



X,,= -t= = 0{\qr-^) (3.23a) 
X.,^-t^ = 0{\qr-') (3.23b) 

All other eigenvalues are Odq]"^"^). 

Let us now return to the case of main interest, in which H and V are the transfer 
matrices (3.2)/(3.3) for the chromatic polynomial (i.e., = — 1 for < i < m, and 

I'j = — 1 for 1 < i < m). In this case we can sharpen (3.15) by providing explicit 
expressions for the lower-order terms: 

Proposition 3.5 Let H and V be the transfer matrices (3.2)/(3.3) for the chromatic 
polynomial = Vi = —1. Then for m> 1, we have that 



t4m) = J2^-l)'ak{m)q'^-' (3.24a) 

k=0 

m 

U{m) = E(-l)V(m)g™-'= + (-1)"^+^ (3.24b) 



fe=0 



^^This Corollary is similar to Corollary 3.3 in [54]. 

^^This Proposition is similar to Propositions 3.4 and 3.6 in [54]. 
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where each akijn) is, for fixed k > 0, the restriction to integers m > max(l,A;) of a 
polynomial in m of degree k given by 



akim 



p=0 ^ 



p(m-p 
P 



k—p 
q=0 



m — 2p 



P 

k — p — q 



(3.25) 



The first few coefficients ak{m) are 

ao(m) = 1 



ai[m) 
a2{m) 
a^im) 



2m + 1 
2m? + m — 2 



fm^ - f m + 1 



(3.26a) 
(3.26b) 
(3.26c) 
(3.26d) 



Tfie proof of Proposition 3.5 is given in Appendix A; it is based on computing tfie 
partition function for a special one- dimensional polymer model. 

Since we have proven that ak{m) is a polynomial in m of degree k, it is also of interest 
to obtain explicit expressions for the coefficients of this polynomial, which we write as 



akim 



{k-e)\{e + 2)\ 



(3.27) 



here the prefactors have been chosen to make many (though not all) of the coefficients 
ak/ integers (in fact, they are all integers for £ < 5 but not for £ = 6, see below). Now we 
use the well-known expansion of the falling powers in terms of Stirling cycle numbers 



X- 



E 

c>0 



(3.28) 



and expand all the binomials in (3.25) involving m. We arrive after some algebra at the 
following expression: 



dki = 



(fc -£)!(£ + 2)! . , . 

- ^ ^ ^ [m" ^]ak{m) 



{k-e)\{i + 2)\{-i) 

2k-2e+i 



(3.29a) 



X 



EE 

a=0 c=0 



^ k k—p 

EE 

p=0 q=0 

k-t 

E 



p 



k—p — qj p\q\ 



-3Y 



d=0 



k-e-d 



r\c—da+c—k+£ 
z p 



(3.29b) 



By computing (3.29b) for integers k > i > 0, we find empirically that a^/ is in fact, for 
each fixed i, (the restriction of) a polynomial in k of degree i. The first few of these 
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polynomials are: 





= 1 




(3.30a) 


«fe,l 


= 3k- 


9 


(3.30b) 


(lk,2 




- 62A; + 4 


(3.30c) 


(lk,3 


= lOA;^ 


- 220A;^ + 410A; + 


(3.30d) 


«fc,4 




- 570k' + 3245fc^ + 1662fc + 2872 


(3.30e) 


«ik,5 


= 21fc^ 


- 1225A;^ + 14245P - Glllk^ + 9982A; + 70560 


(3.30f) 




= 28k^ 


2324A;^ + 45780A;^ 1166396^3 653968^2 








3313360 u 1 1935488 
9 ' 3 


(3.30g) 



Full disclosure. We did not begin by proving Propositions 3.1-3.5. Rather, we 
computed the transfer matrices for small widths m and noticed some patterns. Only then 
did we try to provide proofs. 

3.3 Block-diagonalization of the transfer matrix (II) 

This block-diagonalization can be explained in another way, without reference to trans- 
fer matrices. We begin by observing that the addition-contraction relation (which is a 
variant of the deletion-contraction relation) for multivariate Tutte polynomials gives 

Zsr^Jq,^) = Zs=Jq,v) + Z,^Jq,v) (3.31) 

where 5"^^^ is the graph obtained from Sm,n by contracting the two extra vertices into one, 
while yS"^ „ is the graph obtained from Sm,n by adding av — —1 edge between the two extra 
vertices. Indeed, (3.31) simply says, in coloring language, that the two extra vertices in 
Sm,n must receive either the same color or different colors. The block-diagonalization of T 
into T= and T^ simply corresponds to carrying out separate transfer-matrix calculations 
for the two families 5^ „ and 5*^^ and observing that 

^5=„(g,v) = u;TT=(m)"-ieid (3.32a) 
Zs^Jq,v) = u;^T^(m)--^eid (3.32b) 

Note also that 5*^ „ is equivalent to a lattice of width 171 + 1 with cylindrical boundary 
conditions in which m of the columns have the ordinary vertical edges Vi while one column 
(corresponding to the contracted extra sites) has v — > -|-oo; the horizontal edges are given 
by the as usual (sec Figure 6). Likewise, S^ ,^ is equivalent to a lattice of width 

m + 2 with cylindrical boundary conditions in which m of the columns have the ordinary 
vertical edges Vi while two adjacent columns (corresponding to the two extra sites) have 
V — >■ -|-oo; the horizontal edges are given by the as usual except that the horizontal 
edges between the two adjacent special columns are set to —1 (see Figure 7). 
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For the chromatic polynomials Ps= ^ (q) and P^^ (g), we can use the modified transfer 



matrix T' = HVH because all the horizontal weights are v = —1: 



(3.33a) 
(3.33b) 



The transfer matrix T'^(m) acts on the space of non-crossing non-nearest-neighbor par- 
titions of the set {0, . . . ,m} with periodic boundary conditions (i.e., and m are also 
considered to be nearest neighbors). However, because of the reflection symmetry, we 
can work in the space of equivalence classes modulo reflection (with respect a diameter 
passing through the vertex 0) of non-crossing non-nearest-neighbor partitions of the set 
{0, . . . , m} on a circle. Wc write = dimT'j^(m) for the dimension of this space, where 
T"(m) is the symmetry-reduced transfer matrix for the family Likewise, the transfer 
matrix T^(m) acts on the space of non-crossing non- nearest- neighbor partitions of the set 
{0, . . . , m -|- 1} with periodic boundary conditions. Once again exploiting the reflection 
symmetry, we can work in the space of equivalence classes modulo reflection (with respect 
a diameter passing between the vertices and m+1) of non-crossing non-nearest-neighbor 
partitions of the set {0, . . . ,m + 1} on a circle. We write = dimT'^(m) for the di- 
mension of this space. In Table 1 we quote the dimensionalities N={m) and N^{m) for 
1 < m < 14. 

In Appendix B we compute the following general formulae for A^=(m) and N^{m): 

Theorem 3.6 The number N={m) of equivalence classes, modulo reflection with respect 
an axis going through vertex 0, of non-crossing non-nearest-neighbor partitions of the set 
{0, 1, . . . , m} on a circle is given by 



N={m) 

where m' — [(m -|- 1)/2J and 



lm'/2] 



k=0 



m'\ / m' — k 
k J \k-\- I\m is even] 



(3.34) 



/[condition] 



1 i/ condition is true 
if condition is false 



(3.35) 



The generating function for this sequence is 
G={z) 



1 
z 



-1 

4z 



l-3z 



9- 



1 + z 4z^/{l + z^){l-3z^) 



(3.36) 



In the above theorem, Rjn denotes the mth Riordan number (see Appendix B.l for the 
main properties of the Riordan numbers [9]). We remark that the sequence N={m) coin- 
cides with sequence A005218 in [61] shifted by two units. 
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Theorem 3.7 The number N-^{m) of equivalence classes, modulo reflection with respect 
an axis going between the vertices and m + 1, of non-crossing non-nearest-neighbor 
partitions of the set {0,1, ... ,m,m -\- 1} on a circle is given by 



ik=0 ^ / \ 



m' — k 
-\- I[m is odd] 

where m' — [(m + l)/2j . The generating function for this sequence is 

2z^ - z'^ + 1 



(3.37) 



V(l + ^^)(l-3z2) 



l + z 



2z{l + 2z^) 



(3.38) 



It follows that the dimension of the full symmetry- reduced transfer matrix T"(m) is 

dimT"(m) = SqPree(m + 2) = N^(m) + N={m) . (3.39) 

In particular, by using Theorems 3.6 and 3.7 we have rederived formula (2.17) and have 
found a few equivalent expressions: 

Corollary 3.8 The number dimT"{m) = SqFree(m + 2) of equivalence classes, modulo 
reflection in the center of the strip, of non-crossing non-nearest-neighbor partitions of the 
set {0, 1, . . . , m, m + 1} is given by 



SqFYee(m + 2) = M^+i + - 



2 ^ \ k 

k=0 



k + 1 



-M„« + - [I ) 

i.—n \ / 



k=0 

[m'/2i 



m'\ fm' — k-\-l 
k 



1.. l'V^Um'\fm'-k + 2\ 



(3.40a) 
(3.40b) 
(3.40c) 



where m' — [(m + 1)/2J . The generating function for this sequence is 

^Ii±3iHM-(i+2.+2.=) 



l + z 1 + z"^ 



1 -3z^ 



(3.41) 



We defer to Appendix B the detailed proof of this Corollary and the fact that (3.40a,b,c) 
and (2.17) do in fact coincide. 



3.4 Relating the partition function to the transfer-matrix eigen- 
values 



For each boundary condition ^ (= or the partition function Zgt (q',v) can be 
expanded in the form 



^5i,„(?'V) = ^ak,i{m,q,v)\k,i{m,q,vY , 



(3.42) 
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where Afc^jj(m, g, v) are the eigenvalues (labelled arbitrarily) of the transfer matrix Tn(m, g, v), 
and Q;fc^(|(m, q, v) are the corresponding amplitudes. Please note that the a^,^ and X^^^ are 
algebraic functions of q and v, but we have to be careful about labels as we go around 
branch points. 

For large |g| everything becomes simple. By Corollary 3.4 there exists, for each m 
and V, a value Q{m, v) such that (for both boundary conditions) the dominant eigenvalue 
A^^ji(m,g, v) is analytic in the region \q\ > Q{m,\^) and has there a convergent expansion 
in powers of 

oo 

Kiim, q, v) = ^) ' (3-43) 

fc=0 

with all other eigenvalues being strictly smaller in modulus. We have bl{m,\) = 1 by 
Proposition 3.2 specialized to the particular operators H and V given by (3.2)/(3.3). By 
increasing Q{m,v) if necessary, we may also suppose that the amplitude Q;*,(t(m, g, v) 
corresponding to the dominant eigenvalue is nonvanishing for |g| > Q(m,v). Inserting 
these facts into (3.42), we conclude that for |g| > Q{m,v) we have 

oo 

= l + 5;(-l)^6[(m,v)g-'= 

k=l 

uniformly on compact subsets of |g| > Q{m, v). [Here we used hm„_^oo \V{SLn)\/''^ — "^-l 
It follows from this that (3.44) also holds order-by-order as a power series in 1/q. 

Unfortunately, the foregoing argument does not control the dependence of Q{m,v) 
on m and v. However, the uniformity in m (and v) can be proven, at least for the 
antiferromagnetic case — 1 < I'e < 0, by invoking the following result 

Theorem 3.9 [62, Theorem 6.3] Let G = {V, E) he a loopless finite undirected graph 
in which all vertices have degree < A' except perhaps for an N- clique yi, . . . ,yi\f . Let G 
be equipped with edge weights {fejeeB satisfying —l<Ve<0 for all e and Vl^y^y.) = —1 
for all i ^ j. Then all the zeros of Za{q, v) lie in the disc \q\ < 7.963907 A' -|- N. 

For the family „ we have A' = 4 and = 1, where yi is the extra site; therefore, all 
zeros of Ps=_Jq) lie within the disc |g| < 7.963907A' + 1 = 32.855628. For the family „ 
we have A' = 4 and N — 2, where the 2-chque is formed by the two extra sites; therefore, 
all zeros of (q) lie within the slightly larger disc Igl < 7.963907A' + 2 = 33.855628. 

Using the Beraha-Kahane-Weiss theorem (Theorem 2.2) we can conclude that, outside 
such a disc, the transfer matrix Tjj(m, g,v) for each width m and each antiferromagnetic 
coupling V must have one and only one eigenvalue of largest modulus Av,j(m, q, v). Since 
this dominant eigenvalue cannot collide with any other eigenvalue of the same transfer 
matrix, it must be an analytic function of q (and nonvanishing) outside the given disc; 

-•^^Actually, [62, Theorem 6.3] also gives a more general result applicable to complex edge weights in 
the "complex antiferromagnetic regime" |1 + We| < 1, albeit with slightly weaker bounds. 



(3.44a) 
(3.44b) 
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and the series (3.43) must be convergent outside the given disc. That is, we can take 
Q{m,w) = 33.855628 for all m and all antiferromagnetic v. 

In summary, the transfer matrix for a square-lattice strip of width m and with = or 7^ 
boundary conditions has a single dominant eigenvalue Xi,-{'m,q,v) or X^^^{m,q,v) that 
is a nonvanishing analytic function of q (in fact, g"* times an analytic function of q~^) 
whenever \q\ > 33.855628, uniformly in m and in v belonging to the antiferromagnetic 
regime. 



3.5 Comparing 5^ with 5- „ 

A 2-term,inal graph {G,s,t) is a (loopless and connected) graph G = {V,E) with two 
distinguished vertices s and t (s 7^ t), called the terminals. We define the partial partition 
functions of {G,s,t) as follows [63, Section 2.1]: 

E ^'^^^ (3.45a) 

ACE eeA 
A does not connect s to t 



ACE eeA 
A connects s to t 



Prom (2.4) it follows trivially that 

ZG{q,v) = 4^*H?,v) + 4^*)(g,v). (3.46) 

We define G»st to be the graph in which s and t are contracted into a single vertex. (Note 
that if G contains one or more edges st, then these edges are not deleted, but become 
loops in G • st.) Then from (2.4) it also follows [63] that 

Zcstiq, v) = Zt'\q, v) + q'' zl^^'\q, v) . (3.47) 

Finally, let G + st define the graph G with an extra edge st. Then it also follows from 
(2.4) that 

ZcUq,^) = (l + ^3*)4"*^(?,v) + (l + y) 4^*^(9, v), (3.48) 

which equals Z^q, y)+VstZG»stiq, v) in agreement with the deletion-contraction formula. 

Both ZQ^^\q, v) and ZQ'^*\q, v) are polynomials in q and v. At large the leading 
behavior of Zq"^^^ and Zq'^^^ is 

Z(f^*)(g,v) = + ... (3.49a) 

4-*)(5,v) = E n^^ + --- (3-49b) 

shortest paths eSpath 
s -> t 
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where dcis^t) is the length of the shortest path from s to t using edges having 7^ 0. 

Let us now speciahze these general formulae to the case G = Sm,m with the terminals 
s, t taken to be the two special vertices. Then 



G»st = S, 

G + st ■■ 



m,n 



so that 



(3.50a) 
(3.50b) 

(3.51a) 
(3.51b) 



We have \V{G) \ = mn + 2, \V{G •st) \ — mn + 1 and \ V{G + st)\ = mn + 2. Furthermore, 
we have dcisjt) = m + 1, and there are n shortest paths from s to t: each of them consists 
of the m — 1 horizontal edges in a row of S^^n together with the two additional edges 
joining this row to the vertices s,t. Therefore 



(3.52) 



shortest paths eSpath 
s t 



i=Q 



Finally, in the graph G + st — S^,^ we have Vst--^- 

For \q\ > Q{m,v), using (3.44) applied to S^^ together with (3.51b), we see that 



l/n 



A,,^(m,g,v) = Jim (4^^^*)(g,v)) 

00 

l + 5](-l)\^(m,v)g- 



fc=i 



(3.53a) 
(3.53b) 



Likewise, using (3.44) applied to together with (3.51a), (3.49) and (3.52), we see 
that 



A*=(m,o, v) = hm ( Z^q'^^' 

= A*,^(m,?,v 

= A^,^(m,g,v 

= A*,^(m,g,v 

= A*,^(m,?,v 



(,,v) + ,-4t'(,,v)) 

4r(^>v)i 



lim 

n— J-oo 



1 + g- 



l/n 
l/n 



lim 



1 

m 

l + q-'^nllvi^i+. + Oiq-"^-') 

i=0 

m 

l + g-'"n^MH-i + 0(g-"^-^) 

1=0 

m 

+ l[vi,i+i + 0{q-^) . 



l/n 



(3.54a) 
(3.54b) 

(3.54c) 
(3.54d) 
(3.54e) 



1=0 
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Please note the logic in the step from (3.54c) to (3.54d): we know from (3.44) that the 
limit exists order-by-order as a power series in q~^; so the term of order q~"^ is indeed as 
claimed, and everything else is of order or smaller. We do not need to have any 

explicit bounds on the n-dependence of the 0{q~'^~^) terms in (3.54c). We have therefore 
proven: 

Proposition 3.10 The difference at large \q\ between the leading eigenvalues A* -(m, q, v) 
and A*^^(m, q, v) of the transfer matrices T^(m, q, v) and T'^{m, q, v) for m> 1 satisfies 

m 

A*-(m,g, v) - A*,^(m,g, v) = JJ 0(g"^) . (3.55) 

In particular, for the chromatic-polynomial case Vg = —1 we have 

A.,= (m,g) - A.,^(m,g) = (-1)™+^ + 0{q-') . (3.56) 

The formula (3.55) generalizes the result (2.40) for theta graphs [which corresponds to 
setting the vertical couplings Vi = 0]. It has the same form as the result (3.13) for the 
dominant diagonal entries — the only difference being that (3.13) is exact, while (3.55) 
is for the leading term only. 

4 Numerical results I: Large-g expansion of the lead- 
ing eigenvalue 

In this section we compute the large-g expansion of the leading eigenvalues A* -(m) 
and A*^^(m) and determine empirically some of their remarkable properties. In Section 5 
we shall provide theoretical explanations of some (but not all!) of these empirical obser- 
vations. 

4.1 Overview of method and results 

In Proposition 3.5 (proven in Appendix A) we computed in closed form the dominant 
diagonal entry in the transfer matrix, t= or t^, for a strip of width m > 1 with either = 
or 7^ boundary conditions. We found that this entry is in each case a polynomial in q of 
degree m: 

m 

t^{m) = Y,{-lfak{m)q^-^ (4.1a) 

fc=0 

m 

t={m) = 5^(-l)*^afc(m) g"*-^ + (-1)"^+^ (4.1b) 

We furthermore computed in closed form the coefficients ak{m), which are in fact, for 
each fixed k > 0, the restriction to integers m > max(l, k) of polynomials in m of degree 
k [cf. (3.25)]. 
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In this section we want to carry out an analogous (though ultimately less explicit) 
computation for the leading eigenvalue of the transfer matrix, which we call A^^= or A*_^. 
Prom Section 3.4 we recall that each Aj,(m) has, for large \q\, a convergent expansion in 
powers of q~^, 

oo 

A*,B(m) = Y.{-l)%{m) q-' , (4.2) 

k=0 

where (i denotes = or 7^. It is also illuminating to pass from the eigenvalue Xi,{m) to its 
logarithm, which is a free energy, and define 

^ fe=i 

Corollary 3.4 guarantees that the first two terms in the expansion of the leading 
eigenvalue coincide with those in the dominant diagonal entry: 

b^{m) = b^{m) = ajt(m) for A; = 0, 1 (4.4) 

(with m > 2 required for in order to avoid the (—1)"*+^ term), or more specifically 

6=(m) = 6^(m) = 1 (4.5a) 
b={m) = bf{m) = 2m + l for m > 2 (4.5b) 

Here we shall go further and compute the coefficients h'^{m) and ^^{m) for 1 < m < 12 
and 1 < m < 11, respectively, with < k < AO}^ We shall find empirically that the 
agreement between the eigenvalue and the dominant diagonal entry extends one term 
beyond that guaranteed by Corollary 3.4: namely, we have^^ 

bf{m) — 02(777.) for m > 2 (4.6a) 
6^(m) = a2(m) for m > 3 (4.6b) 

Beyond this, the eigenvalues deviate from the dominant diagonal terms. However, from 
Proposition 3.10 we know — and our computations of course confirm — that the devi- 
ations of the eigenvalues from the dominant diagonal terms are the same for = and 7^ 
boundary conditions, all the way down to order 1, so that the difference between the two 
eigenvalues has the same leading behavior as the difference between the two dominant 
diagonal terms: 

A,,= (m) - A,,^(m) = (-1)™+^ + O(g-i) (4.7) 

as compared with 

t=(m) -t^(m) = (-1)"^+^ . (4.8) 



""^^It would not be difficult to extend this computation to much larger values of k, if we really cared. 
Extension to larger values of m is, however, an extremely demanding computational task. The dimensions 
of the transfer matrices T^(13) and T^(12) are 15465 and 15339, respectively, which are beyond the 
capabilities of our current computer facilities. 

^''In [54] we found, for free and cylindrical boundary conditions, an agreement extending two terms 
beyond that guaranteed by Corollary 3.4. 
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In terms of the coefficients, this says that 



Kim)-bt{m) = \ - - (4.9) 

' -1 tor A; = m 



and 

7^f^\ - Jo forl<A;<m-l 

(-1)'"+^ forA; = m 



c^{m)-ct{m) = \^ (4-10) 



This amazing behavior suggests that we look more closely at the differences between 
= and 7^ boundary conditions for orders k >m. We therefore define coefficients di{m) by 

oo 

(-l)-+MA,,= (m) - A,,^(m)] = 1 + ^(-l)^rf,(m) q-' (4.11) 

or equivalently 

d,{m) = - [6=+,(m) - C(m)] . (4.12) 

It furthermore turns out to be interesting to look at the logarithm of this difference of 
eigenvalues: 

oo 

log{(-l)'"+^[A.,= (m) - A.,^(m)]} = ^e,(m) q-' . (4.13) 

e=i 

The following variant is also of interest: 



log 



(-l)-+^^-log(^^^)l = f2fi(m)q-r (4.14) 



i=i 



We shall obtain empirical formulae for the coefficients di{m), e^(m) and /^(m) for 1 < 

e < 13. 

In this section we shall proceed as follows (imitating what we did in [54, Section 4] for 
free and cylindrical boundary conditions): First we shall compute the transfer matrices 
for strips of width m < 11 (resp. m < 12) for 7^ (resp. =) boundary conditions, using 
the methods of [34,53] as modified for the family Sjn,n (see Section 3 above). Prom these 
matrices we then extract the dominant eigenvalue as a power series in q~^, i.e. for each 
available m we compute as many coefficients b^j^{m) and c[,(m) as we please. -"^^ We then 
observe empirically that, for each A; > 0, the coefficients b'^{m) and bf{m) [resp. c^{m) 
and c^(m)] are equal and are represented by a polynomial Bk{m) [resp. Cfc(m)] of degree 
k [resp. degree 1] in m, provided that we restrict to integers m > mmin(A;), where 

^min(^) — max(/c, 1) (4.15a) 

^min(^) ^ k + 1 (4.15b) 



^^To compute the dominant eigenvalue as a power series in q~^, we have used the power method [28, 
Section 7.3.1] in symbolic form. Each iteration gives one additional term in the expansion of the dominant 
eigenvalue in powers oi q~^. We can therefore compute the exact expansion up to any desired order in a 
finite number of steps. 
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for the two boundary conditions. ^'^ Assuming that this empirical observation is accurate 
(i.e., that the polynomial behavior observed for m < 11, 12 persists to all larger m), we 
can infer the exact expressions for the polynomials and for /c < 17 for 7^ and = 
boundary conditions. The fact that the polynomials Cfc are of degree 1 in m — hence 
consist of a bulk term akm and a surface term Pk — can be explained (non-rigorously) 
by finite-size-scaling theory: see Section 5 below and especially the discussion around 
eqns. (5. 7)/ (5. 9) and at the beginning of Section 5.5. 

Analogously — and rather more surprisingly — we find that the coefficients rf^(m), 
e^(m) and fi{m) are also represented by polynomials D^{m), E^{m) and F^{m), of degrees 
^, 1 and 1, respectively, provided that we restrict to integers m > I + 1. This behavior 
— in which the difi^erence of two eigenvalues (which are essentially free energies) behaves 
like a partition function in the sense that its logarithm grows asymptotically linearly in 
the strip length — is wholly unexpected, at least to us. 

4.2 Results for ^ 

Using the methods just described, we have obtained the leading eigenvalue A*^=^(m) 
for 1 < m < 11 as a power series in [cf. (3.43)] through order k — 40. The resulting 
coefficients bf{m) are displayed for < A; < 15 in Table 2, and the corresponding coeffi- 
cients c^{m) [cf. (4.3)] are displayed in Table 3. (The complete data set for < A; < 40 
is contained in the Mathematica file data_Neq.m that is included in the on-line version 
of this paper at arXiv.org.) It is interesting to note that for all pairs {k,m) that we have 
computed (namely, 1 < m < 11 and < k < 40), the coefficients b^{m) and kc^{m) 
are integers. Furthermore, we observe empirically that, for each fixed k, the coefficients 
6^(m) [resp. c^(m)] are the restriction to integers m of a polynomial Bk{m) [resp. Cfe(m)] 
of degree k (resp. degree 1) in m, provided that we restrict attention to m> rri^^^{k) with 

"^min(^) = rnax(/c, 1) . (4.16) 

Below this threshold ni^^^^^^k), the coefficients deviate from polynomial behavior. With 
our available data together with a few tricks described below, we are able to determine 
these polynomials for < k < 17. 

First we start by trying to fit the coefficients b'^{m) with m > m^^^{k) to a polynomial 
Bk{m) of degree k. As we need k + 1 coefficients for such a polynomial, we are able to 
obtain these polynomials when rn^^^^{k) + A; < 11, i.e. A; < 5. Please observe that in all 
cases we have at least one data point more than the number of unknowns, so that every 



^''By contrast, for the dominant diagonal entry we have proven in Proposition 3.5 that ak{m) are 
polynomials in m of degree k; and in this case the polynomial form holds for all allowable integers m, 
i.e. TO > max(fc, 1). 
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fit can be tested on at least one extra data point. Our results are: 



Boim 

Bi{m 
B2{m 
Bs(m 



= 1 



2m + 1 
2m^ + m — 2 
- f m + 3 



- Im^ - + - 3 



15 



(4.17a) 
(4.17b) 
(4.17c) 
(4.17d) 
(4.17e) 
(4.17f) 



We see that the three highest-order coefficients in these polynomials agree with those of 
the corresponding polynomial ak{m) [cf. (3.25)/(3.26)], i.e. 



ak{m) 



forO<A;<2 

Bk{m) = { 

ak{m) + 0{'nr ^) for /c > 3 



(4.18) 



There is, however, a better way of extracting the desired information from our numer- 
ical data: instead of using the coefficients ^^{m) as our basic quantities, we can use the 
related coefficients c^{m) [cf. (4.3)]. The latter coefficients are empirically found to be, 
for each fixed k, the restriction to integer m of a polynomial Ck{m) of degree 1 , i.e. 

Ckirri) = akUi + Pk , (4.19) 

provided again that we restrict attention to m > the same m^^^i^k) defined in (4.16). As 
we now need only the two coefficients and (3k to specify such a polynomial, we are able 
to obtain these polynomials up to A; = 9 (if we want at least one extra data point to test 
the fit) or A; = 10 (if we don't). The first polynomials are given by 



Ci(m) = 


-2m - 1 


(4.20a) 


C2(m) = 


-m-| 


(4.20b) 


C3(m) = 


>-f 


(4.20c) 


C^m) = 


5 41 
2"i 4 


(4.20d) 


C5(m) = 




(4.20e) 



and the results for 6 < Ac < 10 are shown in Table 4. The polynomials for A; < 10 can 
then be determined from the Ck using (3.43)/ (4.3). 



^®In our view this test is important, because the formula (4.16) for rn^i^{k) is only empirical; and if 
it were to be mistaken by one unit (i.e. the real threshold for polynomial behavior were one unit higher 

than wc thought), then a fit with no extra data points would give a definite result (one can always fit 
k + 1 points to a polynomial of degree k) but this result would be total nonsense. 
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Actually, we can do better than this. We believe (or at least conjecture!) that the 
coefficients c^{m) are, for each fixed /c > 0, the restriction to integers m > m^^^{k) of a 
polynomial Ck{m) of degree 1. If we now compute the difference 

Atim) = ct{m)-C,{m) (4.21) 

between the numerical coefficients d^{m) and the corresponding polynomials Ck{m), we 
find, not surprisingly, that these differences are nonzero whenever m < m^i^(/c): see 
Table 5, where we initially know A^(m) only for A; < 10 (or < 9 if we wish to be 
conservative). If we could somehow guess an analytic form for at least some of these 
coefficients A^(m), we could then define "improved" coefficients c^{m) by 

cf(m) = 4{m)-At{m), (4.22) 

so that these coefficients ^{m) would be, for each fixed k, the restriction to integers m > 
^min(^) '^^ ^he same polynomial Ck{m), but with a smaller threshold rn^^^(k) < iTi^-^^^{k). 
The important point here is that a smaller threshold ^min(^) imphes that we can obtain 
more polynomials Ck with the same raw data. 

By inspecting Table 5, it is not difficult to realize that there are some patterns in 
A^(m) just below the threshold rri^^^{k), and more specifically along diagonals m — k — £ 
with fixed integer i > 1. For instance, ior i — 1 we see that A^(A; — 1) = (—1)*^+^ for 
2 < A; < 10; and it is reasonable to conjecture that this holds also for larger k, i.e. 

A^(A;-1) = (-1)^=+^ forallA;>2. (4.23) 

Likewise, for £ = 2, 3 we find for A; < 10 — and conjecture for all larger k — that 
A^(A; — i) can be written as (—1)'^+^ times a polynomial in k of degree i — 1, provided 
that k>2e-l: 

Af(A;-2) = (-1)'=+^ (A; - 7) for A; > 3 (4.24a) 

A^(A;-3) = (-1)^+^ (iyt' - f A; + 57) for A; > 5 (4.24b) 

Each of these fits can be tested on at least three additional data points. [We can also 
obtain a fit for A^(A; - 4) for A; > 7, but at this stage we would have no extra data points 
to test it.] 

Then, assuming the correctness of (4.23)/(4.24), we can infer the polynomials Ck for 
k — 11, 12 (if we demand at least one extra data point to test the fit) and A; = 13 (if we 
don't). In turn, knowing the polynomials Cn and C12, we can obtain a fit for A^(A; — 4) 
that can be tested on two extra data points: 

Af{k-i) = - we + ^k - m) for A; > 7. (4.25) 

Next, assuming the correctness also of (4.25), we get one additional data point to test 
the fit to Ci3 (not surprisingly, it checks) and we can infer C14 if we don't demand any 
extra data points to test it. In turn, knowing the polynomial C13, we can guess one further 
correction term A^(/c — £): 

Atik - 5) = (-1)^+^ (^A;^ - §k' + - + 4253) for > 9 , (4.26) 
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where here there are no additional data points to test the fit if we use only C13, but 
there is one additional data point if we trust also the inferred C14. Finally, assuming the 
correctness of (4.26), we can infer C15 (but without any extra data points to test the fit). 
All the polynomials for A; < 15 are shown in Table 4. 

We can extend slightly our results by noticing the following two empirical results from 
the formulae for A^(/i; — C) in (4.23)-(4.26): (a) the leading coefficient in A^(A; — £) is 
l/(£ - 1)!, and (b) the next-to-leading coefficient is -l/(£ - 1)! x [0, 7, 27, 60, 106], which 
is fitted by |(^ — 1)(13^ — 12). If these patterns persist, then we need only 4 coefficients 
(rather than 6) to get A^(/c — 6), i.e. k — 11, 12, 13, 14, so we can then get C15 with one 
extra data point, or Cie without. The result of the fit is 

(A;-6) = (-1)^^+1 (^fc^ - ^A;^ + iff - ^k'' + »A; - 378651) for > 11 , 

(4.27) 

where here there is one additional data point if we trust also the inferred C14 and C15. 

Remark. In Section 5.5 we will infer also C17 by exploiting the known large-g series 
for the bulk free energy [1,54], which gives us the coefficients a^; we then need only one 
additional data point to infer 13k and hence Ck- AH the polynomials Ck for A; < 17 are 
shown in Table 4. 

4.3 Results for ^ 

We carried out a completely analogous analysis for the case of = boundary conditions; 
we shall therefore be brief in describing it, stressing only the differences from the preceding 
analysis and the final results. 

We were able to compute \^-{m) for for 1 < m < 12 (one term more than for 7^ 
boundary conditions). The resulting coefficients b^{m) and c^irn) for < A; < 15 are 
displayed in Tables 6 and 7, respectively.^^ We observe empirically that, for each fixed 
k, the coefficients b'^{m) [resp. c^{m)] are the restriction to integers m of a polynomial 
Bk{m) [resp. Cfe(m)] of degree k (resp. degree 1) in m, provided that we restrict attention 
to m > m^^^{k) with 

^min(^) ^ k + 1. (4.28) 

Furthermore, the polynomials Bk{m) and Ck{m) are empirically found to be the same as 
those for 7^ boundary conditions. Using our available data together with the tricks shown 
in the preceding subsection, we are able to confirm this equality for A; < 17. Let us sketch 
quickly the logic of this computation and the relevant intermediate results. 

As before, we can extract the desired information most efficiently by fitting the coef- 
ficients c^(m) to a polynomial Cfe(m) of degree 1, for m > m^^^{k) = A; -|- 1. As we need 
only two coefficients for such a polynomial, we can obtain these polynomials up to A; = 9 
(if we want at least one extra data point to test the fit) or A; = 10 (if we don't). The 
results are equal to those of 7^ boundary conditions. 

Next we try to improve these results, as before, by defining the difference A^(m) = 
c^(m) — Cfc(m) and attempting to guess an analytic form for some of the coefficients 

^^The complete data set for < /c < 40 is contained in the Mathematica file data_Eq.iii that is 
included in the on-line version of this paper at arXiv.org. 
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A^(m): see Table 8, where we initially know only k < 10. As in the case of 7^ bound- 
ary conditions, we find empirically that the coefficients A^(m) closest to the boundary 
"^min(^) restriction to integers m of certain polynomials. Because m^^^{k) = k + 1, 

we are concerned this time with A^(/c — i) ior i > (rather than ^ > 1 as before). We 
obtain: 



Ar(fc) = 


(- 


-1)^+^ for > 1 


(4.29a) 


AH^-i) = 


(- 


-1)^+^ {k - 5) for A; > 3 


(4.29b) 




(- 


.l)fe+i(iA;2_MA; + 37) for A; > 5 


(4.29c) 


Ar(fc-3) = 


(- 


-1)'=+^ (Ik^ - 9^ + ^k- 303) for A; > 7 


(4.29d) 



We are able to test each fit on at least one additional data point. The pattern is 
slightly different from that for 7^ boundary conditions: here A^(A; — i) is (—1)'^"'"^ times 
a polynomial of degree £ in k, vahd for A; > 2£ + 1. By this method, we can obtain the 
polynomials up to k — 13 (if we demand at least one extra data point to test the fit) 
or k — 14: (if we don't). 

We can then guess one further correction term A^(A; — i): 

A^(A;-4) = (-l)'=+i (^A;^ - f|Ar^ + ^A;2 - i^A; + 2573) for A; > 9 . (4.30) 

where there arc no additional data points to test the fit if we use only C13, but there is one 
additional data point if we trust also the inferred C14. Finally, assuming the correctness 
of (4.30), we can infer C15 (but without any extra data points to test the fit). 

As for ^ boundary conditions, we can extend slightly our results by looking at the 
patterns for the first two terms in A^(A; — £) (4.29d)/(4.30): The leading coefficient in 
A^(A; — f.) is l/f.\, and the next-to-leading coefficient is —l/£\ x [0,5,23,54,98], which is 
fitted by ^£{13£ — 3). If these patterns persist, then we need only 4 coefficients (rather 
than 6) to get A^(A; — 5), so we can then get C15 with one extra data point, or Ciq without. 
The result of the fit is 

A^(A;-5) = (-1)^-^ (-^P + ifc^ - i|f fc^ + '-^k' - + 22330) for > 11 , 

(4.31) 

where here there is one additional data point if we trust also the inferred C14 and C15. 
We find, in all cases, the same polynomials Ck as for 7^ boundary conditions. 

Remark. In Section 5.5 we will infer also C17 for = boundary conditions by exploiting 
the known large-g series for the bulk free energy [1 , 54] . 

4.4 Comparing S^^^ with S=^^ 

We were able to compute the difference (— l)™'+^[A^^=(m)— A+^^(m)] for for 1 < m < 11, 
and we have extracted the coefficients d^{m) [cf. (4.11)] for 1 < m < 11 and 1 < £ < 20. 



^•^For f = 3 we do not immediately have such an extra data point, but we can obtain it after using 
using £ = 0, 1, 2 to infer Cn and C12. The reasoning is completely analogous to what was explained in 
detail for the case of ^ boundary conditions. 
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The resulting coefficients are displayed in Table 9 for 1 < £ < 13. (The complete data set 
for 1 < £ < 20 is included in the Mathematica file data_Dif f .m that is included in the 
on-line version of this paper at arXiv.org.) We observe empirically that, for each fixed i, 
the coefficients d(,{m) are the restriction to integers m of a polynomial D^{m) of degree I 
in m, provided that we restrict attention to m> mmin(£) = £ + 1 [cf. (4.28)]. We start by 
trying to perform such fit: As we need i+1 coefficients for the polynomial Df {m), we are 
able to obtain these polynomials for £ < 4 (if we allow for at least one extra data point 



to test the fit) or i <5 (if we don't). Our results are: 

Di{m) = 3m -4 (4.32a) 

D2(m) = Im'^-fm + U (4.32b) 

D3{m) = - f + 84m - 55 (4.32c) 

D,{m) = fm^-lfm3 + lflm2-ifim + 224 (4.32d) 

D^im) = §m5 - i|im^ + ^m^ - ^m^ + l^m - 880 (4.32e) 



As in the preceding sections, we can extract more information by considering the 
coefficients e^(m) arising in the logarithm of the difference of eigenvalues [cf. (4.13)]. 
These coefficients are displayed in Table 10forl<£< 13, and the complete set is 
included in the file data_Dif f .m. We again observe empirically that, for each fixed i, the 
coefficients e^{m) arc the restriction to integers m of a polynomial Ef irn) of degree 1 in 
m, provided that we restrict attention to m > mmin(£) = i+l- As we need two coefficients 
for such a polynomial, we can obtain these polynomials up to ^ = 8 (if we want at least 
one extra data point to test the fit) or £ — 9 (if we don't). The first polynomials are given 

by 



Ei{m) = 


-3m + 4 


(4.33a) 


E^im) = 


-^m + 6 


(4.33b) 


E^{m) = 


-16m + f 


(4.33c) 


Ei{m) = 


-ifim + 66 


(4.33d) 


E^{m) = 


658™ I 934 
5 5 


(4.33e) 



and the results for 6 < £ < 9 are shown in Table 4. The polynomials D£ for £ < 9 can 
then be determined from the Eg using (4.11)/(4.13). 

Next we try to improve these results, as before, by defining the difference 

(m) = ee{m)-Ee{m) (4.34) 

and attempting to guess an analytic form for some of the coefficients Af\k): see Table 11, 
where we initially know only £ <9. As in the previous cases, we find empirically that the 
coefficients Af\k) closest to the boundary A;min(-^) — £ + 1 are the restriction to integers £ 
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of certain polynomials. We obtain: 

Af(£) = (-1)^ for£>l 



-1 



-9) for £ > 3 



Af(£-2) = {-ly {-e^ - 8e + 39) iori>5 



Af\e-3) = i-iy {ff - 65f + - Q13) ior£>7 



(4.35a) 
(4.35b) 
(4.35c) 
(4.35d) 



We are able to test each fit on at least one additional data point. By this method, we are 
able to extract the polynomials up to £ = 12 (if we demand at least one extra data 
point to test the fit) or £ = 13 (if we don't). 

The same game can be played with the coefficients fe{m) [cf. (4.14)]. These coefficients 
are displayed in Table 12 for 1 < ^ < 12, and the complete set is included in the file 
data_Dif f .m. From this table we see that the coefficients lfi{m) are not integers, in 
contrast to the observed behavior for the coefficients lei{m) (see Table 10). We observe 
empirically that, for each fixed £, the coefficients /^(m) are the restriction to integers m of 
a polynomial F^{m) of degree 1 in m, provided that we restrict attention to m> rrimmi^) — 
£+1. Again, we can obtain these polynomials up to £ = 8 (if we want at least one extra 
data point to test the fit) or £ = 9 (if we don't). The first polynomials are given by 



Fi{m) 


~ m — 5 


(4.36a) 


F2{m) 


= >-¥ 


(4.36b) 


Fs{m) 


_ 49_ 77 
~ 3 3 


(4.36c) 


F4(m) 


_ 191 _ 305 
— 4 4 


(4.36d) 


F5{m) 


= ^m-203 

5 


(4.36e) 



and the ones for 6 < £ < 9 are shown in Table 4. (In Section 5.4 we prove that fi{m) 
m — 5 for m > 2.) 

We now try to improve these results by defining the difference 



A'^/\m) = Mm)-Fe{m) 



(4.37) 



and attempting to guess an analytic form for some of the coefficients A^^-* (m) : see Table 13, 
where we initially know only £ <9. As in the previous cases, we find empirically that the 
coefficients A^/\m) closest to the boundary minin(-^) — £+1 are the restriction to integers 
£ of certain polynomials. We obtain: 

; > 1 (4.38a) 
L) for £ > 3 (4.38b) 
! - 54) for £ > 5 (4.38c) 
^-f +^-fi- 736) for £ > 7 (4.38d) 



A'-/\£) = 


(-1)^+^1 f. 


Aif\£-1) = 




a'^/\£-2) = 


i-iY {\e + 


A''/\i-3) = 


i-iy (f £3 - 
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We are able to test each fit on at least one additional data point. By this method, we are 
able to extract the polynomials up to £ = 12 (if we demand at least one extra data 
point to test the fit) or i — 13 (if we don't). 

If we look carefully at Table 4, we reahze that the polynomials and F( are not 

independent, but they satisfy (at least for 1 < £ < 13) the relation 



Ci{m) = Ee{m) + Fi{m). 



(4.39) 



The same relation of course holds for the corresponding coefficients cf{m), e^(m) and 
fe{'m) for m > £ + 1: 

cf{m) = ee{m) + U{m) for m > £ + 1 . (4.40) 
We can understand this latter relation by using (4.3)/(4.13)/(4.14) and defining 

^\ 

e^im) + fi{m) - cf{m) 



G{m) = 



e=i 



log 



/ K-{m 
X 



A*,= (m) 



(4.41) 



Since A* -(m)/A^,^(m) = 1 + 0{q-'^) from (3.43)/(4.5a)/(4.7), it follows that G{m) = 
0(5"™), i.e. \q~^\G{m) = for £ < m, which is precisely (4.40). 



4.5 Summary of conjectured behavior 

Let us conclude this section by summarizing our empirical findings concerning the 
eigenvalues A*-(m) and A*^^(m) and their difference. First recall that in Proposition 3.10 
we prove that 

A,,= (m) - A,,^(m) = (-1)-+^ + 0{q-^) . (4.42) 

We have found empirically that stronger results hold Our first finding concerns the be- 
havior of the free energies for ^ and = boundary conditions: 

Conjecture 4.1 There exist polynomials Ci, C2, . . . with rational coefficients, of degree 1, 
such that 

log^^F^ = ^C,(m)g-^ + 0(g-(-+^)) (4.43a) 
^ fe=i 

log = ^Cfc(m)g-'= + (-l)"^+^g-™ + ©(g-^'^+i)) (4.43b) 

^ fc=i 

We have verified this conjecture for m < 10, using the polynomials Ci, . . . , Cio shown in 
Table 4. A special case of Conjecture 4.1 is Proposition 3.10 concerning the difference of 
eigenvalues. 

Secondly, the empirical observation that the coefficients di{m) defined in (4.11) are 
the restriction to integers m > £ + 1 of a polynomial of degree I in the variable m can 
formalized as a conjecture as follows: 
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Conjecture 4.2 There exist polynomials Di, D2, . . . with rational coefficients, with deg D£ — 
i, such that 



A. 



■1) 



m+l 



m— 1 



-lYDe{m)q- 



(4.44) 



We have verified this conjecture for m < 6, using the polynomials -Di, . . . , given in 
(4.32). Moreover, the same polynomials Di, . . . , give correctly the expansion through 
order q"^ for all m < 11. 

Finally, the empirical observation that the coefficients etirn) and fiim) defined in 
(4.13)/(4.14) are the restriction to integers m > ^ + 1 of a polynomial 0/ degree 1 in the 
variable m can formalized as a conjecture as follows: 



Conjecture 4.3 There exist polynomials Ei, E2, . . . and Fi, F2, . . . with rational coeffi- 
cients, of degree 1, such that 

m—l 

log [(-l)-+MA,,=(m) - A,,^(m))] = J] £;,(m)?-^ + 0(?— ) (4.45) 

e=i 

m— 1 

= 5^F,(m)g-^ + 0(g-'") (4.46) 
e=i 

We have verified this conjecture for m < 10, using the polynomials Ei,...,Eg and 
Fi, . . . , Fg shown in Table 4. 



-log 



:-i 



im+l m 



g^log 



5 Thermodynamic limit m ^ cx) of the free energies 

In the previous section we studied the largc-g expansion of the leading eigenvalues 
Xi,-{m) and A^^^(m) — or what is essentially equivalent, the strip free energies and 

to be defined in (5.2) / (5.4) below — for strips of fixed width m with = or 7^ boundary 
conditions. In this section we will study the thermodynamic limit m — )■ cxd of these free 
energies. 

The plan of this section is as follows: In Section 5.1 we introduce some preliminary 
definitions and discuss the expected behavior of the strip free energies per unit length, 
fmi^) fm{'i)i ^ function of the strip width m. We then give, in Section 5.2, an 
overview of the computations we have performed and the results we have obtained. In 
Section 5.3 we present a hand computation of the first four terms in the large-q' expansion 
of the bulk, surface and corner free energies for = and ^ boundary conditions. In Sec- 
tion 5.4 we prove some exact results, valid to all orders in 1/g, for the difference between 
the free energies with = and ^ boundary conditions. Finally, in Section 5.5 we explain 
one of the empirical observations from Sections 4.2 and 4.3, and we extract additional 
information on the vertical-surface free energies from our transfer matrices. 
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5.1 Generalities and finite-size-scaling theory 

Corollary 3.4 shows that, for each width m and each boundary condition (= or 7^), 
the transfer matrix has, for sufficiently large \q\, a single dominant eigenvalue A^(g) that 
moreover is an analytic function of q (in fact, it is q"^ times an analytic function of q~^). 
Theorem 3.9 actually proves that X-k,]j,{q) is q"^ times an analytic function of q~^ whenever 
\q\ > 33.855628, uniformly in m. 

Let us now introduce the free energy per site^^ for a finite strip with = or 7^ boundary 
conditions, 

where ji denotes = or 7^, and its limiting value for a semi-infinite strip, 

fjq) = lim — logP^i (g) . (5.2) 

Finally, let us introduce the free energy per site for the infinite lattice, 

f{q) = lim — logP^s {q) . (5.3) 

Here we are assuming that the indicated limits exist and that in (5.3) the limit is inde- 
pendent of the way that m and n tend to infinity. Furthermore, it is natural to expect 
that in (5.3) the limiting free energy is independent of boundary conditions, in which case 
we can omit the superscripts = or 7^ and write simply f{q). 

In fact, some of these assumptions can be proven. Indeed, the above discussion guar- 
antees that, at least for \q\ > 33.855628, the hmiting strip free energy fm{q) exists for all 
m and is given by 

= ^^^SKM, (5-4) 

fit 

which is moreover an analytic function of q in the indicated domain. Moreover, Procacci 
et al. [49, Theorem 2] have proven that, when \q\ is large enough (namely, l^l > 8e^ ^ 
160.684295), the infinite-volume limiting free energy f{q) exists and is analytic in 1/g and 
is the same for all sequences of graphs Gmxn with free boundary conditions; in particular, 
it is independent of the way that m and n tend to infinity. We expect that the same 
result holds for = or 7^ boundary conditions and gives rise to the same infinite-volume 
limiting free energy /(g) as for free boundary conditions (indeed, we expect that this can 



Actually minus the free energy per site in the usual thermodynamic convention — but we prefer not 
to encumber our formulae with unnecessary minus signs, or to encumber our text with constant repetition 
of the word "minus" . 

^^In fact, Procacci et al. prove that this limit is the same for all F0lner van Hove sequences of finite 
subvolumes of the infinite lattice, not just rectangles; and they prove this result for a wide variety of 
lattices (namely, locally-finite connected quasi-transitive amenable infinite graphs), not just the square 
lattice. On the other hand, we mistakenly asserted in [54, Section 5.1] that Procacci et al. had proven 
this result also for cylindrical, cyclic and toroidal boundary conditions. We suspect that such a result 
can indeed be obtained by a suitable modification of their proof, but we were wrong to assert that it is 
contained in their paper. 
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be proven by a modification of the Procacci et al. argument). In this paper we will take 
n ^ oo first and then take m ^ oo, so that 

m = lim fl{q) . (5.5) 

Finite-size-scaling theory [48, Section 2.5] gives a rather precise prediction for the form 
of the free energy (5.1)/(5.2) for a finite or semi-infinite system away from a critical point 
(and in the absence of soft modes). In particular, for an m x n strip with = or 7^ boundary 
conditions and bulk correlation length ^buik <^ m,n, the predicted behavior is 

71 f'^ ~\~ Til f'^ 1 

ftt — £ „ _L surf, vert ^ ./surf.horiz _|_ .[1 , (^( -m.m{m,n)/^bulk\ (r. a) 

Jm,n — ./bulk -r mn ^ ^ ' ^ 

where ^ denotes = or ^, and /bulk = /, if, vert' /sLf,horiz and /J^.^,, arc, respectively, the 
free energies for the bulk, the two vertical surfaces (at left and right), the two horizontal 
surfaces (at top and bottom), and the four corners plus the extra sites. For a semi- infinite 
strip m X oo with = or 7^ boundary conditions, we have 

fm = /bulk + ;^/Lf,vert + ©(e""^/^--) . (5.7) 

lit 

Since the families S^^^ are not invariant under 90° rotation, we expect the vertical and 
horizontal free energies to be different. Finally, we cannot take for granted that the surface 
and corner free energies are the same for = and 7^ boundary conditions. 

In this section we assume that the behaviors (5.6)/(5.7) hold and that expansion in 
1/m and 1/n can be commuted freely with expansion in 1/q. (We expect that this may 
be provable by an extension of the Procacci et al. argument.) This assumption justifies 
the manipulations to be made in the following subsections. 

The relation (5.7) of course holds for the chromatic polynomials at fixed large q. But 
we can also argue heuristically what it should imply for the series expansion in powers of 
1/q. It is not difficult to see that, for large g, we have 

e-i/«buik(g) ^ I + o(^) (5.8) 
q \q''/ 

(just as for a one- dimensional Potts antiferromagnet at zero temperature). As explained 
in detail in [54], we can therefore interpret 0(e~"*/'^''""=^*^) as meaning 0{q~^). Therefore, 
we expect that 

fUq) = /bu.k(g) + ^/Lf,vert(?) + 0{q-n , (5.9) 

or in other words we predict ni\^:^^{k) = k + 1. Of course, we should not take too seriously 
the here, since the amplitude of the correction term in (5.7) could be proportional 

to a positive or negative power of q. But we do predict that ni\^:^JJz) kin the sense that 
limfe_^oomJjji^(A;)/A; = 1. This is indeed what we found for = and 7^ boundary conditions. 
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5.2 Overview of computations and results 

Before beginning the detailed computations, it is useful to give an overview of the 
methods to be used and the results obtained, taking into account both the computations 
reported below for = and 7^ boundary conditions and those reported in a previous paper 
[54, Sections 5.2-5.4] for free (F) and cylindrical (P) boundary conditions. 

The computations we are able to carry out are as follows: 

• For the bulk free energy (which is of course the same for all four boundary condi- 
tions), we obtained in [54] long series (through order g^"^^) using the finite- lattice 
method for free boundary conditions. We are also able to check these series up to 
moderate order (namely, order q~^^, q~^^, q~^^ respectively for F, P, = and 7^ 
boundary conditions) by transfer-matrix calculations of the leading eigenvalue. 

• For the vertical surface free energy with = or 7^ boundary conditions, we have a 
hand calculation through order q~'^, which we can extend to order q~^^ or q"^^ by 
transfer-matrix calculations of the leading eigenvalue. 

• For the horizontal surface free energy with = or 7^ boundary conditions, we have 
only a hand calculation through order g^^. [The horizontal surface free energy does 
not appear in the strip eigenvalue (5.4)/(5.7).] 

• For the comer free energy with = or 7^ boundary conditions, we have likewise only 
a hand calculation through order q~^. 

Our conclusions (from the available orders) are as follows: 

• The bulk free energy agrees for all four boundary conditions (no surprise). 

• The horizontal surface free energy agrees for = and 7^ boundary conditions. This 
is again no surprise, as the horizontal surfaces are the same for the two boundary 
conditions. 

• The vertical surface free energy agrees for = and 7^ boundary conditions. A priori 
this might be somewhat of a surprise, as it is far from obvious whether the effect 
of the edge connecting the two extra sites (which is w = cxd or — 1 according as the 
boundary conditions are = or 7^) is of order 1 (it is, after all, only a single edge) or is 
of order n (its endpoints are directly connected to the whole vertical sides) . It turns 
out that the former is the case. But on closer examination this is no surprise at all: 
the equality /subvert — /surf vert follows immediately from (5.4)/(5.9) together with 
Proposition 3.10, since log[A*-(m)/A^,,5^(m)] = ( — l)™+^g^™ -|-0(g^™^^) vanishes to 
all orders in q~^ as m — >■ 00. 

• The horizontal and vertical surface free energies are different from each other, and 
different from the surface free energy for free (F) boundary conditions. 

• The corner free energies are different for =, 7^ and free boundary conditions. How- 
ever, the difference between the corner free energies for = and 7^ boundary conditions 
is simple: f^mevil) - ftomeM = -log(l - q-^). 
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Let us conclude this subsection by recalling the large-g series for the bulk free energy 
that was obtained in [54] through order g""^^ (extending previous series by Bakaev and 
Kabanovich [1]) by using the finite-lattice method [22-25,43,46] for free boundary con- 
ditions.^^ As |g| — > oo, the exponential of the free energy per site for an infinite square 
lattice is given by the series expansion 

^ (izil! [1 + ^ + ^ + 3^8 + 4/ + 3^10 + 3^11 + 11^12 + 24^13^8^^^ 
Q 

-91z^^ - 261z^^ - 290^^^ + . . . -598931311074^^^ + 0{z^^)] , (5.10) 
where z is defined as 

(5.11) 

We have here copied the first few terms and the last one; the remaining terms are reported 
in [54, Table 9]. In terms of the variable 1/q, we obtain 

e-^(9) ^ q[l- 2q-^ + q'^ + q'^ + q'"^ + q'^ + q'^ + 2q-^ + 9q-^ + SSq'^ + 130?"^° 
+378g-^^ + 987g-^2 + . . . + 1311159363081366872?-^^ + 0{q-^^)] . (5.12) 

Finally, for the bulk free energy f{q) itself (rather than its exponential) in terms of the 
variable 1/q, we obtain 

1 I 1 I 5 I 28 I 55 I 89 I 81 I 505 I 1029 , 7742 



f(n\ — Intr ,7 _ 2 _ 1 I 1 I 5 I 28 I 55 I 89 I 81 I 505 I 1029 , 7742 

I 25291 I 73552 , 197755 , , 190018276619486037135 , /0^„-48^ /'cioN 

1"T2^ + T3?3 + + ■ ■ ■ iV'T \-^{Q )■ \0.L6) 



Clearly, e-^*^'^-' has a much simpler expansion than f{q)] in particular, its coefficients are 
integers (at least through the order calculated thus far).^^ A further simplification is 
obtained by using the variable z — l/{q — 1) in place of 1/q: the integer coefficients 
become much smaller. Finally, a slight extra simplification arises from extracting the 
prefactor (g — l)^/g in (5.10). 

Unfortunately, the finite-lattice method — which is an extraordinarily efficient method 
for calculating series expansions — does not appear to be applicable to = and ^ boundary 
conditions. We therefore resorted to a hand calculation, to be reported in the next 
subsection, which we carried through order q~'^. 

5.3 Hand calculation of large-g expansion for the bulk, surface 
and corner free energies 

First of all, as we are interested in the large- limit, it is convenient to explicitly 
remove the leading term logg' in the free energy by considering the modified chromatic 



After the completion of this work, we learned that Jacobsen [32] had extended the series expansions 
for the bulk, surface and corner free energies to orders 0(g~^^), 0{q~^^) and 0{q~'^^), respectively. 

^*The coefficients of f{q) are not integers, but k[q~'^]f{q) is an integer (at least through the order 
calculated thus far). Indeed, it is not hard to show that if F{z) is a power series with integer coefficients 
and constant term 1, then k [z*^] logF(2:) is always an integer. 
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polynomial Pq for a loopless graph G — {V,E): 



(5.14) 



Using the Fortuin-Kasteleyn representation (2.4) we get 



^kiA)-\V\ 



(5.15a) 



ACE 



l)|A|^l/^)|A|-c(A) 



(5.15b) 



ACE 



where c{A) = \A\ — \ V\ + k{A) is the cyclomatic number of the subgraph (V, A). 

It is instructive to begin by computing "by hand" the first few terms of the large-gr 
expansion for the bulk, surface and corner free energies. To do this, let us first consider 
an m X n square lattice with ^ boundary conditions: it has \V\ = mn + 2 sites, \E\ = 
2mn — m + n+l edges, |Fsq| = (m — l)(n — 1) square faces, and |Ftri| = 2(n — 1) triangular 
faces. There are also two extra faces, each bounded by a cycle of length m + 2, cf. the 
inner and outer faces in Figure 7b; but these faces will play no role in the calculation if 
m is large enough (see below). We can compute the first few first terms in the large-g 
expansion for the modified chromatic polynomial P^j^ (q) by using (5.15b) and explicitly 
identifying the subsets A having a given small value of \A\ — c{A): 



\A\ - c{A) = 0: Only A^0. 

\A\ — c{A) ~ 1: A — any single edge. 

\A\ —c{A) =2: A — two distinct edges or three edges forming a triangular 
face. 

1^41 — c{A) — A — three distinct edges not forming a triangular face, 
or four edges forming a cycle of length four (i.e. a square face or the 
boundary cycle of two adjacent triangular faces), or five edges consisting 
of all the edges of two adjacent triangular faces. 



We have also identified the contributions of the subsets with \A\ —c{A) = 4 and will include 
them in our calculation, but we refrain from writing them down as their description is 
rather lengthy. For small m we also have to worry about terms A that wind horizontally 
around the lattice using the two extra sites; the smallest of these are cycles of length 
m + 2. But since all such terms have l^l — c{A) > m + 1, we can avoid them simply 
by assuming that m is large enough, i.e. m > k if we want an expansion valid through 
order q~''. In particular, we can obtain the expansion through order by assuming that 
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m > 4 and ignoring "winding" subsets A. We therefore have 



1^1(1^1 -1) 



-Pt.ri 



1 



+ 



\E\{\E\-V){\E\-2) 
6 

1^1(1^1 -1)(|^|-2)(|E|- 3) 



|Fsq|-|Ftri|(|i?|-2) 



24 



(|FtH| + |F,q|)(|£;|-3) 



and hence 



+0(g 



— l-^tri 

5\ 



(|^|-3)(|E|-4) , iFtriKlFtril-l) 



+ 



(5.16) 



1 

2^ 



2 J^l - 2|Ftri 



1 

3g' 



3 jE|-3|F,q|+6|Ftri 



4g 



^(|£;|-12|Fsq| + 14|Ftri|) 



(5.17) 



Note that in the expression for logP^^^ (q'), all geometrical quantities occur linearly (as 
they should). Inserting the values of \E\, \F^^ and \F^^-\, dividing by mn, expanding in 
1/n and 1/m, and putting back the leading term logg, one finds the large- expansion for 
the free energy 



/m,n(?) 



2 115 
logg ^ + + + 



1 

9 



5 



16 

3^ 



41 



1 

m 



+ 



1 1 
g 2g2 

+o(g-^). 



2 

3^ 



11 



+ 



1 3 14 39 



1 

ran 



Comparing to the finite-size-scaling Ansatz (5.6), we obtain 

/buik(g) 



2 1 1 5^,5, 



^ 3urf,vert(^) 



^ surf jhoriz (^) 



5 



^6 
3^ 



4g^ 



(5.18) 

(5.19a) 
(5.19b) 
(5.19c) 
(5.19d) 

mn + 1 sites, 

(m — l)(n — 1) square faces, and |Fti.i| = lin — 1) 
triangular faces. There are also two extra faces, each bounded by a cycle of length m+1, cf. 
the inner and outer faces in Figure 6b; but these faces will play no role in the computation 



f corner l^) 



15 16 41 ^, _5, 

q^ 2q^ 3g3 4g4 + ^ 
1 3 14 39 ^, 5, 



We can do the same computation for the family S^^: it has \y\ 



\E\ = 2mn — m + n edges, \F, 



sq I 
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if m is large enough. The subsets A contributing through order 0(g~^) are exactly the 
same as for the family S^^. The only difference is that the "winding" terms now start 
with cycles of length m + 1 rather than m + 2, and they have \A\ — c{A) > m instead 
of m + 1. Therefore, we must take m > A; + 1 if we want an expansion valid through 
order q~^. In particular, we can obtain the expansion through order by assuming that 
m > 5. The large-g expansion for logPs= Jg) is given by the same formula (5.17) but 
with the new value of inserted (note that |Fsq| and \Firi\ remain the same). Dividing 
by mn, expanding in 1/n and 1/m, and putting back the leading term log q, one finds the 
large- expansion for the free energy 



fm,niQ) = logg---- + 

1 1 2 

q 



1 5 
3^^ 



1 



5 



16 
3^ 



41 



m 



11 



1 


' 2 


5 


10" 


1 


- + 










n 


g2 






mn 



+0(g-^). 

Comparing to the finite-size-scaling Ansatz (5.6), we obtain 

logg 



/bulk(g) 
3urf jvert (^) 
f 3urf,horiz 
f comer (?) 



2 11 5 ^, 



1 

2 



9 

5 

2^ 
1 



16 
2 

3^ 



41 



4g^ 
11 



+ 0(g-^) 



4g' 



,4 + 0(?~') 



5 10 _5. 



(5.20) 

(5.21a) 
(5.21b) 
(5.21c) 
(5.21d) 



If we compare the series (5.19) with (5.21), we conclude that: 

1) The bulk contributions are the same for ^ and = boundary conditions, and they 
agree also with the result for free and cyhndrical boundary conditions [54, Eq. (5.14a)]. 

2) The surface contributions are the same for both boundary conditions, and they de- 
pend on the orientation of the "surface" (i.e. vertical or horizontal). Moreover, /surf honz ~ 
ffati horiz ^-g^^^es with the surface free energy for free boundary conditions [54, Eq. (5.14b)]. 
This is not unexpected, as the top and bottom rows are identical for all these boundary 
conditions. 

3) The corner free energies do depend on the boundary conditions (i.e. /comer 7^ 
/corner) ^^^l they also differ from the corner free energy for free boundary conditions [54, 
Eq. (5.14c)]. We observe the curious fact 

/c=omer(?) " /ctner(?) = ^ + ^ + ^ + + ^(^"') " (5-22) 

In the next subsection we will establish this equality to all orders in 1/g. 



5.4 Exact results for the bulk, surface and corner free energies 

We can obtain some exact results concerning the difference between the free energies 
fmrSs) ^^"^ frnrSs) — ^^"^ heucc between the bulk, surface and corner free energies for 
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= and 7^ boundary conditions — by using the results of Section 3.5. Let us write the 
modified chromatic polynomials Pot (q) [cf. (5.14)1 in terms of the partial chromatic 
polynomials of the two-terminal graphs {Sm,n, s,t) introduced in Section 3.5; we obtain 



-,—mn—l 



,-1 pisi^t). 



q 



Ojn..n. N -*■/ 



(5.23a) 
(5.23b) 



Therefore, the difference between the finite-size free energies for the two boundary condi- 
tions can be written as 



fmAq) - ftniq) 



mn 



^ogPs= (q) - log P^^iq) 



mn 



log 1 - - H log 



mn 



1 + q 



P, 



(q) 



(5.24a) 
(5.24b) 



where P^*^*-* and P^*^*-* 



are the corresponding modified chromatic polynomials. From 
(3.49) and (3.52) we see immediately that 



(q) 



1 + 0{q-') 



n 



and hence 

fZ,n{q) - /m,n(?) 



mn 



log 1 - 



+ 



\m+l 



m 



(5.25a) 
(5.25b) 

(5.26) 



The term proportional to q ^ disappears to all orders inl/gasm— >oo, hence we have 



f 3urf ,horiz(?) 
f surf, vert (5) 

f corner (?) 



f 3urf,horiz(?) 
f surf,vert(?) 



/corner (?) " log (^l - 



(5.27a) 
(5.27b) 

(5.27c) 



to all orders in 1/q?^ In particular, this proves (5.22) to all orders in 1/q. 

We can actually push the expansions (5.25a,b) to higher order in 1/q by methods 
similar to those used in the preceding subsection. Indeed, Ps\^^^ is exactly equal to P^ ^ 
whenever m is large enough (compared to the order in 1/q being considered) so that 
winding subsets A (i.e. sets connecting the two terminals s, t) are absent. In particular, 
PsS'n^ is given through order q~'' {k < 4) by (5.16) whenever m > k. On the other hand, 
a straightforward hand computation yields 



P, 



(?) 



:-i) 



m+l 



7m+l 



n — 



n\E\ + mn — 5n — m + 4 



+ 0{q-') 



(5.28) 



25 



We can interpret the correction term 0{q as 0(e '"/^buH') [cf. (5.8)]. 
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1 / p^'""'\Qy 



whenever m > 2. It follows that 

log A* -(m) - log A*,^(m) = hm log | 1 + g ^^^^^^ 



n->oo n 



q" 



m — 5 



(5.29a) 



(5.29b) 



In terms of the quantities fe{m) defined in (4.14), this says that fi{m) = m — 5 for m > 2, 
which agrees with Table 12. 

5.5 Families S^ ,^ and S^ ,^ via strip free energies 

Using finite-size-scaling theory, we can understand (non-rigorously) a fact observed 
empirically in Sections 4.2 and 4.3: namely, that the coefficients c^(w^) arising in the large- 
\q\ expansion of the free energy (4.3) are represented for m > fnl^^^{k) by a polynomial of 
degree 1 in the strip width m, i.e. 



Cfe ("^) = cf(m) = Ck{m) = akin + Pk for m > mmin(A;) . 



(5.30) 



To see this, it suffices to compare the large- jg] expansion of the hmiting free energy for a 
semi-infinite strip. 



fmiQ) = -logA.,tt(g) 
m 

^ oo 

= log g+— log ^{-l)''bl{m)q 
Lfc=o 

1 °° 

- logg+- J](-l)'=4(m)g-\ 



(5.31) 



fc=i 



with the finite-size-scaling Ansatz (5.7)/(5.9). The behavior (5.30) is an immediate con- 
sequence, where (resp. (5k) is the coefficient of q^^ in the bulk (resp. vertical surface) 
free energy. As discussed after (5.9), we expect ml^^Jyk) fa k. 

We can also use our transfer matrices to check in part the result (5.13) for the bulk 
free energy and to notably extend the results (5.19b)/(5.21b) for the vertical surface free 
energy. Indeed, the coefficients in both these expansions can be immediately read off from 
the coefficients ak and (3k in Table 4. We thus obtain an independent check of the ffist 
15 terms of the series (5.13). We also obtain the first 15 terms in the large-g expansion 
of the vertical surface free energy /^tf.vert = /surf , vert ^ 



surf,vert 



1 

<1 



49 55 719 2459 

3^ V 8^ "9^ 
1239 15168 23051 171677 647719 744743 



5 



16 



41 



81 

bq^ 



2gio llgii 
6898415 ^, 



6g 



12 



13g 



13 



Uq 



14 



bq 



15 



16g 



16 



(5.32) 
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We can slightly extend this latter series by using the series (5.13) for /bulk as an input: 
in this case, each polynomial Ck contains a single unknown coefficient to be determined 
(rather than two unknown coefficients). We then obtain the coefficient of the term 
in f'^ 

J surf, vert' 

[Q-''] /irf,vert = ^ ' (5-33) 

and thus the polynomial C17 shown in Table 4. 

Finally, the simplest expression for the surface free energy is given by its exponential 
in terms of the variable z — 1 / {q — 1): 

e4Vf,vert(9) = 1 _ 2 _ ^2 ^ ^4 ^ 3/ + 2/ - 5^^ - 11^^ + 3/ + 43z^° + 57 z'^^ - SAz^^ 
-178^13 - 122^^^ + 220^^^ + 200^^*^ - 1170^^^ + 0{z^^) . (5.34) 



6 Numerical results II: Limiting curves Bm 

For each 1 < m < 10, we have computed the symbolic transfer matrix T"(m), as 
well as the vectors fjd and a;""". In this computation we exploited the fact that T"[m) is 
block- diagonal [cf. (3.31)/(3.33)], which allows us to obtain separately the two diagonal 
blocks T^(m) and T^(m) and, for each block, the corresponding left and right vectors 
oj"^ and fid. 

The symbolic computations were done by means of a Mathematica program for 
1 < m < 5 and by a perl script for 1 < m < 10. The perl script runs faster and uses a 
smaller amount of memory than the Mathematica program for the same width m. Wc 
have checked that both programs give the same answer for 1 < m < 5. We also performed 
several checks on our symbolic results using the identity 

Psr^jQ) = Ps^Jq) ' (6.1) 

where Sm,n is the graph that is just like Sm,n except that the two extra sites lie at top and 
bottom rather than left and right (see Figure 8). Of course Sjn,n and Sn^^ are isomorphic, 
so that the identity (6.1) holds trivially, but the transfer- matrix approaches to calculating 
their chromatic polynomials are rather different, because of a reversal in which direction 
is considered "longitudinal" and which "transverse" (see Appendix C). The identity (6.1) 
thus constitutes a nontrivial check on the correctness of our computations. We checked 
it for 1 < m < 8 and 1 < n < 6. 

For each 1 < m < 6, we computed the zeros of the chromatic polynomials Psm.pmi^) 
for strips of aspect ratio p = 10, 20, and wc also computed the accumulation set Bm of 
chromatic roots in the limit p — )■ 00. For m < 4, wc used the resultant method [53, Section 
4.1.1] to compute the limiting curve Bm, together with the direct-search method [53, 
Section 4.1.2] to refine some details. For m = 5, we used the resultant method to compute 
the endpoints of the limiting curve, while the rest of the curve was obtained via the direct- 
search method. For m — 6, the full limiting curve was obtained using the direct-search 
method. We took advantage of a few slight improvements to these methods as described 
in detail in Ref. [38, Section 2]. In principle we could have determined the limiting curve 
Bm also for 7 < m < 10; but it is extremely laborious (in both human and CPU time) 
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to determine accurately the small details in Bm using the direct-search method, so we 
decided to stop at m = 6. 

Since H is a projection, the basis elements in our transfer- matrix calculation are given 
by fp = Hep, and not by e-p, as explained in Section 2.2. However, to lighten the notation, 
wc will represent the basis for a given strip as a collection of vectors e-p; but it should 
be understood that the basis vectors are actually f-p. We shall use the delta-function 
shorthand (2.13) to write the vectors e-p. 



6.1 m = 1 

The transfer matrix T'^(l) is one- dimensional, and a basis is given by B= — {!}. This 
matrix and the corresponding left and right vectors are, in this basis, 

r^(l) = q-2 (6.2a) 
uj"^ = q{q-l) (6.2b) 
= 1 (6-2c) 

The transfer matrix T'^(l) is also one-dimensional, and a basis is given by = {I}- 
In this basis, the transfer matrix and the corresponding vectors are given by 

7^(1) = q-3 (6.3a) 
ojT^ = g(g_l)(g_2) (6.3b) 
= 1 (6.3c) 

The two eigenvalues are evidently X= — q — 2 and = ? — 3, which become equimod- 
ular along the line 

Reg = ^. (6.4) 

Thus, this is the limiting curve for this strip. The amplitudes can be easily read off from 
(6.2)/(6.3): i.e., a= — q{q — 1) and a^t — q{q — l){q — 2), in agreement with the exact 
result (2.28). 

The dominant eigenvalue in the half-plane Reg < 5/2 is = g— 3, while the dominant 
eigenvalue in the half-plane Reg > 5/2 is A= = q — 2. The amplitude of the dominant 
eigenvalue vanishes at g = 0, 1, 2 (and only there); hence, these three points are the only 
isolated limiting points for this strip. The limiting curve Bi and the chromatic roots for 
n = 10, 20 are depicted in Figure 9. 



6.2 m = 2 

The transfer matrix T'^(2) is one-dimensional. In the basis B= — {1}, this matrix 
and the corresponding vectors take the form 

T=(2) = g2-5g + 7 (6.5a) 
ojT ^ ^(^_i)(^_2) (6.5b) 

f5 = 1 (6.5c) 
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The transfer matrix T'^(2) is two-dimensional. Let us choose the basis as = 
{^0,2 + ^1,3, !}• The transfer matrix and the left and right vectors for this strip are then 

u;T _ (2(g-l),g2_3g + 3) (6.6b) 

= (0,1) (6.6c) 
The eigenvalues are given by 

A= = g^-5g + 7 (6.7a) 

A^,i = l(l0-6q + q^-y^^ (6.7b) 




A^,2 = ^ (lO - 6g + g2 + ^g^ j (6.7c) 
where we have used the shorthand notation 

Q2(q) ^ q^- 8?^ + 28?' - 56g + 52 . (6.8) 
The eigenvalues A^,i and A^,2 are the solutions of the quadratic equation 

- x{q^ - 6g + 10) - + Iq^ - 16g + 12 = . (6.9) 
The amplitudes are given by 

= g(g_l)(g_2) (6.10a) 
q{q - 1) (3 - 3g + q^)^2 - (22 - 34g + 21g' - + g^) 

2 yg^ ^ ■ 

g(g - 1) (3 - 3g + q^)^2 + (22 - 34g + 21g' - Iq^ + g^) 

"^'^ ^ ^ 7m ^ 

We have computed the limiting curve B2 by using the resultant method. This curve 
crosses the real g-axis at g = 3, which is a quadruple point. There are two pairs of complex- 
conjugate T points, namely q ^ 2.6099757836 ± 1.8423725343 i and q ^ 2.8900242164 ± 
0.5194968788 i. 

In the two regions having nonempty intersection with the real g-axis, the dominant 
eigenvalue comes from T'^(2); in the other four regions, the dominant eigenvalue is A=. On 
the complex-conjugate curves connecting the T points q 2.6099757836 ± 1.8423725343 i 
and q ^ 2.8900242164 ± 0.5194968788 i, A^,i and A^,2 are equimodular. Hence, at the 
four T points the three eigenvalues become equimodular. Finally, at the quadruple point 
A= = A^,2 = 1 and A^,i = 0.'^ 



^^The eigenvalues A^^i and A^,2 are analytic functions of q except at the branch cuts of the function 

1/(32(9) • The polynomial Q2 has zeros at q = 2 ± i + =F 2i « 2.78615 T 0.27202i and at g = 
2=Fi — V— 1 ± 2i 1.21305=F 2.27202i. These four zeros belong to regions where A= is dominant: namely, 
the first pair of zeros belong to the two oval-like closed regions protruding from the quadruple point g = 3, 
while the later pair belong to the other two regions not intersecting the real q-axis. We have chosen the 
branch cuts to be horizontal lines starting at each zero and going to g — >■ — 00: i.e., for each zero qi with 
i = 1, . . . , 4, we define q = qi+ riC^^* with > and — tt < 6i < tt. With this definition of branch cuts, 
A^^2 is dominant in the two regions having nonempty intersection with the real g-axis, except whenever 
2.27202 > I Img| > 0.27202 and Req lies to the left of the lines joining two T points. 
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There are four isolated limiting points at g = 0, 1, 2, B^. At these values the dominant 
amplitude q;^^2 vanishes. The limiting curve B2 and the chromatic roots for n = 20, 40 
are depicted in Figure 10. 



6.3 m = 3 



The transfer matrix T'j^(3) has dimension 3. The basis is chosen as B= = {5o,2, ^1,3, !}• 
In this basis, the transfer matrix and the left and right vectors for this strip are 





/ -L4A -1 


T:(3) 


= { q-2 












= (0,0,1) 



5-2q 
1 



where we have used the shorthand notation 

Ti{q) 



- Iq^ + 19g - 20 



M) 



q — mq + n 



(6.11a) 

(6.11b) 
(6.11c) 



(6.12a) 
(6.12b) 



The transfer matrix T'_^(3) has dimension 4. The basis is chosen as — {^0,2 + 
^4,2, (^0,3 + (^1,4, (^1,3, 1}- In this basis, the transfer matrix and the left and right vectors for 
this strip are 



r'(3) 



id 



/ -L5.6 g-3 -1 3-g\ 

q-2 -L5,7 q~2 1 

q-2 1 

\ 2L5 6 2L5 7 Lg^g Ti — 1 J 

= qiq -l)iq- 2) {2{q - 1), 2(g - 1), g - 1, L,,^) 
= (0,0,0,1) 



(6.13a) 

(6.13b) 
(6.13c) 



where we have used (6.12). 

The three eigenvalues of T^(3) come from the third-order polynomial 



X 



- x^{q^ - 8g2 + 24g - 26) - x{q^ - I2q^ + 59?^ - 149?^ + 193g - 101) 



- IZq^ + 70?^ - 2mq^ + 320?^ - 212q + 96 (6.14) 
while the four eigenvalues of T^(3) come from the fourth-order polynomial 

- a;^(g^ - 9g2 + 30g - 36) -x2(2g^-26g^ + 140g''^- 388^2 + 551g- 318) 



X 



-x{q'^ - 19q^ + 153g^ - 681g^ + 1815g3 - 2901^^ + 2577g - 981) 



+q^ - 19g^ 



157g^ - 738g^ + 2161g^ - m9q^ 



+4707g2 - 3128g + 908 . 



(6.15) 



We have computed the limiting curve by using the resultant method, and fine-tuned 
the results using the direct-search method. The limiting curve B crosses the real q- 
axis at g ~ 2.8177131118. There is one pair of complex-conjugate endpoints at g ~ 
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2. 7852013976±0. 7798713630 i. There are seven pairs of complex- conjugate T points at g ^ 
1.0229779807 ± 2.4615628248 i, q ^ 2.3154363287 ± 1.9575822608 i, q ^ 2.4538824566 ± 
1.8334262878 q ^ 2.5862741033 ± 1.5583461255 i, q ^ 2.8988477958 ± 0.4358504561 i, 
q ^ 2.9168615472 ± 0.6195940190 i, and q ^ 3.0476880109 ± 0.4560842031 i. There are 
four isolated limiting points at q — 0, 1, 2, B5. 

The dominant eigenvalue comes from T'^(3) in the region containing the interval q G 
(—00,2.8177131118), and in the two large regions bounded asymptotically by |argg| e 
(7r/6,7r/2) as \q\ — )• 00. The dominant eigenvalue comes from T'^(3) in the other three 
large regions: the one that contains the interval q e (2.8177131118, 00) and the other two 
asymptotically bounded by | argg| e (7r/2, 57r/6) as \q\ -^00. 

There are four isolated limiting points at g = 0, l,2,i?5. The limiting curve and 
the chromatic roots for n = 30, 60 are depicted in Figure 11. 

6.4 m = 4 

The transfer matrix T"(4) has dimension 5, and T'^(4) has dimension 9. They 
are too lengthy to be reported here, but they can be found in the Mathematica file 
square_extra_sites.m (1.9 MB) that is available on request from the authors. 

We have computed the limiting curve by using the resultant method, and fine-tuned 
the results using the direct-search method. The limiting curve B crosses the real q-axis 
at q ^ 2.9060325277 and q ^ 3.9030119682. There are two pairs of complex-conjugate 
endpoints at g ^ 2.6169563471 ± 1.4357301392 i, and q ^ 2.8849105593 ± 0.7356893954 i. 

There are ten pairs of complex-conjugate T points atq ^ 0.3411402631±2. 3118521368 i, 
q ^ 1.1260056618±2.5364550333i, q ^ 1.3027336503±2.5447465991 i, q ^ 2.2640047808± 
1.9520686076 i, q ^ 2.4920548627 ± 1.7925248200 i. q ^ 2.7002452362 ± 1.5923412229 i, 
q ^ 2.8652031850±1.153161817H,g 2.9601783483±0.5518584736 i, g 2.9623803372± 
0.5807622346 i, and q ^ 3.0624715624±0.6434188421 i. There is a pair of bulb-like regions 
protruding from T points q ^ 2.7002452362 ± 1.59234122291 

The dominant eigenvalue comes from T^^(4) in the region containing the interval q e 
(—00, 2. 9060325277)U(3. 9030119682, 00), and in the two large regions bounded asymptot- 
ically by I arggi G (37r/8, 57r/8) as |g| — )> 00. The dominant eigenvalue comes from T'j^(4) 
in the other four large regions: i.e., those asymptotically bounded by | argg| G (7r/8, 37r/8) 
and I argg| G (57r/8, 77r/8) as |g| — )■ 00. The dominant eigenvalue comes from this block 
also in the region containing the interval q G (2.9060325277, 3.9030119682), and the other 
two closed regions pointing to the right in Figure 12. 

There arc four real isolated limiting points at g = 0, l,2,i?5, and a pair of complex- 
conjugate isolated limiting points at g ~ 2.8555521103 ± 0.9018551071 i. The limiting 
curve B4 and the chromatic roots for n — 40, 80 are depicted in Figure 12. 

6.5 m = 6 

The transfer matrix T'^(5) has dimension 11, and the matrix T'^(5) has dimension 21. 

This is the first case where we have computed the limiting curve using the direct-search 
method, except for the endpoints that were computed using the resultant method. The 
limiting curve B crosses the real g-axis at q 2.9268787368. 
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There are 14 pairs of complex-conjugate T points ai q ^ 0.0265421747±2. 1428032766 i, 
q fti 0.2211937394±2.2964381892^, q ^ 0.5678532323±2.4981722822 i, q ^ 1.5459127222± 
2.4751057958 q ^ 1.5503506823 ± 2.4742675615 i, q !^ 2.2165122888 ± 2.0475480936 i, 
q ^ 2.3467238502±1.9855375152i, q ^ 2.5384138871±1.7717517039i, q ^ 2.8417052689± 
1.2459939540 i q ^ 2.9278218101 ± 0.9693176128 i, q ^ 2.9651600535 ± 0.7574018277i, 
q ^ 3.0055390406±0.9222249388i,g ^ 3.1196012894±0.266671987H, and g 3.1587639645± 
0.5586982282 i. There are two pairs of complex-conjugate endpoints at g ~ 2.0144392334± 
2.0277231065 i, and q ^ 2.3246979483 ± 1.90163769711 Finally, there is a pair of 
complex-conjugate bulk-hke regions protruding from the T points q ^ 2.3467238502 ± 
1.98553751521 

There are four real isolated limiting points at g = 0, 1, 2, i?5, and two pairs of complex- 
conjugate isolated limiting points at q f=i 2.7190757419±1.4455587779 i, and q a; 2.8265910048± 
0.9420673312 i. The hmiting curve B5 and the chromatic roots for n — 50, 100 are depicted 
in Figure 13. 

6.6 m = 6 

The transfer matrix T'j^(6) has dimension 21, and the matrix T'^(6) has dimension 49. 

The limiting curve has been completely determined by using the direct-search method. 
This curve B crosses the real g-axis at g a 2.9477131589, and g a 4.2138764783. 

There arc 16 pairs of complex-conjugate T points at g ~ 0.1195227877±2.2972440367i, 
g fti 0.8945659676±2.5558662562«,g 1.0187320850±2.5519732169 g fa 1.7361257994± 
2.4200684634 g ^ 2.0496846729 ± 2.2109675624 i, g « 2.0659399006 ± 2.1928587989 i, 
g ^ 2.1816875305±2.1274251306i, g ^ 2.5037573375±1.8512451004i, g ^ 2.7553322546± 
1.4607519460 i, g ^ 2.9072366988 ± 1.0781535615 i, g ^ 2.9302627055 ± 1.0265712687i, 
g ^ 2.9319763271±1.016109882H,g ^ 2.9600871035±0.8562446926 g ^ 3.0829802406± 
0.7932348764 i, q ^ 3.0859998597±0.4641385413 i, and g ft; 3.0934600966±0.4747067803 i. 
There are two pairs of complex-conjugate endpoints at g fs 1. 8355583156±2. 1088360158 i, 
and g fa 2.8978250409 ± 1.0665199267i. 

In this case we have been able to find four isolated limiting points at g = 0, 1, 2, B^; 
but we cannot rule out the existence of other isolated limiting points (especially complex 
ones). The limiting curve Be and the chromatic roots for n = 60,120 are depicted in 
Figure 14. 

6.7 Comparison of 1 < m < 6 

In Figure 15 we show the limiting curves Bm for 1 < m < 6, plotted together. We see 
clearly that: 

(a) There is an oval-shapcd region at small g, extending roughly from g = to g = 3 in 
the real direction and between g 1.1 ± 2.6 i in the imaginary direction, where the 
limiting curves Bm do not enter, or at least from which they retreat as m ^ 00. 

(b) Outside of this region, the curves Bm appear to become dense as m ^ 00, except 
possibly near the real axis. 
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(c) The curve Bm also exhibits, for each m > 3, m — l small "fingers" extending outside 
the oval-shaped region towards the right (these "fingers" are more easily seen on 
Figures 11-14). 

In Figure 15 we have also superposed, in dark gray, the limiting curve for a 
square-lattice strip of width m = 11 and cylindrical boundary conditions [34]. This curve 
can be taken as a rough approximation of the expected infinite-volume limiting curve 
'Soo(sq) = lim for the square-lattice model with cylindrical boundary conditions, 

m— >-oo 

with the exception that we expect the curves B^' to move slightly outward as m grows 
and to close up at g = 0; in particular we expect Soo(sq) to be a closed curve that crosses 
the real g-axis at g = and go = 9c = 3 [34,53]. We expect, furthermore, that the same 
limiting curve will be obtained for free boundary conditions [34, 53] . 

Figure 15 clearly suggests that the oval-shaped curve to which the curves Bm are 
retreating as m — >■ oo is the same as the curve Boo(sq) obtained for cylindrical (or free) 
boundary conditions. The region where the curves Bm are apparently becoming dense is 
precisely the exterior of B^x) (sq) . 



7 Conjectures on the limiting curves Bm as m — oo 

In this section we discuss in more detail the behavior at large |g| of the limiting curves 
Bm', in particular, we formulate conjectures concerning the behavior of these limiting 
curves as m — > oo. 



7.1 Behavior at large \q\ for each m 

Figures 9-14 show the hmiting curves Bm for 1 < m < 6: in each case we give both 
a "close-up" view that shows the details of the small- |g| behavior and a "distance" view 
that makes clear the large-|g| asymptotics. We notice, first of all, that the curve Bm has 
2m outward branches tending to g = oo, with asymptotic (as |g| — ?> oo) angles that are 
equally spaced around the circle. These empirical findings for 1 < m < 6 motivate the 
following result: 



Proposition 7.1 Fix m > 1. Then the limiting curve Bm for chromatic roots of the 
square- lattice strips Sm,n (n ^ oo) has exactly 2m outward branches extending to infinity, 
with asymptotic angles arg g = 6k,m where 

Ok,m^{k-]^- fork = l,...,2m. (7.1) 
\ 2 J m 

Moreover, the dominant eigenvalue comes from TL{m) in the asymptotic regions 

. /Iff A; = 1. 3. . . . , 2m — 1 if m is even ^\ 

^.,^<argg<^,+,,^ for | ^ ^ 3, 4, . . . , 2m if m is odd ^^'^^ 

while in the other asymptotic sectors the dominant eigenvalue comes from T'!^ (m) . 
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Proof. By the Beraha-Kahane- Weiss theorem (Theorem 2.2), the hmiting curve B 
consists of those q for which 

A^-(m) 



log 



A^-(m) 



1, 



or equivalently 
or equivalently 

L ^KA^) 
[the factor ( — is included here solely for future convenience]. If we now define 

A*-(m) 



-l)-'log(^) 



i9 with 9 real , 



TT 

(2A; — 1)77 with /c integer 



(7.3) 
(7.4) 

(7.5) 



^m(g) = -log 



(-l)"*+^?"*log 



[compare (4.14)], then (7.5) can be rewritten as 



TT 



KA"^) 



(7.6) 



argg = (2A; — 1) Im 

2m m 

On the other hand. Proposition 3.10 together with (4.5a) tells us that 

•^m(g) = 0{q-^) , 

so that the curve Bm is given by 

argg = (2fc-l)^ + 0{\q\-^) . 

This proves the first statement of Proposition 7.1. 

The second statement of Proposition 7.1 follows easily from 



(7.7) 



(7.8) 



(7.9) 



which implies 



KA^) 



l)™+^g-™ + 0(g-("^+i)) 



KA^) 



= 1 + (-l)"*+'Re(g-"^) + 0(g-("^+^)) 



(7.10) 



(7.11) 



□ 



Remcirk. This proof of Proposition 7.1 actually uses a bit less than the full strength 
of Proposition 3.10. More specifically, the first statement of Proposition 7.1 follows from 



A*-(m) — A^,^(m) = Cm + o(l) with Cm real 7^ 



(7.12) 



while the second statement uses sgn(c„i) = (—1)™+-'^. In neither case do we need the exact 
value Cm = (—1)™^^, nor do we need to know that the error o(l) is actually 0{q~^). □ 
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Using Proposition 3.10 we can actually go farther, and compute the corrections to 
argg = 9k,m as a series in inverse powers of It suffices to use the fundamental relation 
(7.7) and insert the expansion 



(7.13) 



[cf. (4.14)], which is guaranteed by Corollary 3.4 together with Proposition 3.10 to be 

convergent for all sufficiently large \q\ (how large may depend on m). It is then straight- 
forward to determine iteratively argg as a power series in For instance, the series 
through order \q\~^ is given by 



argg 



e 



k,m 



+ 



+ 



fiim) 



sin 6, 



m 

\3 



k,m 



+ 



h{mf , /2(m) 
2m? 



+ 



m 



sin 29k \ 



fiimy , /i(m)/2(m) 



2m? 



(3 sin 39k,, 



sm 9k,m) H sm36'fe,^ 



m 



+o{\qn. 



(7.14) 



In particular, for m < 11 we know explicitly the coefficients fi{m). 

If we assume Conjecture 4.3, then we can replace /^(m) by the polynomial F^{m) 
whenever m > i + 1, and we can moreover use the specific forms of F^{m) shown in 
Table 4 at least for £ < 9. 

In Figures 9-14 we compare the exact limiting curves Bm with the approximations 
(7.14) through orders and q~^. We have checked numerically the correctness of 

the asymptotic expansions (7.14), by plotting the deviations of &xgq from (7.14) truncated 
at order \q\^^ ■ niTiltipIicd by versus and verifying that they are tending to 

constants (or at least bounded) as \q\ oo. 



7.2 Uniformity in m and density of chromatic roots 

Since the asymptotic rays argg = 9k,m arising in Proposition 7.1 become dense as 
m — 7> oo, it is natural to expect that the chromatic roots of the graphs Sm,n become 
dense as m, n ^ oo in the whole complex plane outside of a bounded set (analogously to 
what happens for the generalized theta graphs [63]). We can formalize this conjecture as 
follows: 

Conjecture 7.2 There exists a constant Q < oo such that the chromatic roots of the 
graphs Sm,n become dense in the region \q\ > Q when m,n ^ oo. 

We stress that Conjecture 7.2 docs not follow from Proposition 7.1 alone — or from the 
more fundamental Proposition 3.10 and Conjectures 4.2-4.3 — because these results apply 
to each value of m separately, and there is no guarantee that the error bounds in these 
asymptotic expansions are uniform in m; in particular, there is no guarantee that they 
can be made to apply to a region 1^1 > Q where Q is independent of m. In this subsection 
we would like to discuss a refinement of Proposition 3.10 from which Conjecture 7.2 can 
indeed be deduced. 
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We know rigorously that there exists Qq < oo, independent of m, such that A*^=(m) 
and A*^^(m) are analytic and nonvanishing in \q\ > Qq for all m > 1 (see Theorem 3.9 
and the discussion following it, where it is shown that Qo — 33.855628 suffices). If we 
now use Proposition 3.10, we can conclude that the functions J-mio) defined in (7.6) are 
likewise analytic in |g| > Qo for all m > 1 and moreover vanishing at g = oo, so that they 
are given for |g| > Qo by a convergent expansion (7.13). Let us now make the following 
hypothesis: 

Conjecture 7.3 There exist Qi e {Qq, oo) and C < oo such that 

sup |J>n(q')| < Cm for all m . (7-15) 

l9l=Ql 

Note that here Qi and C are independent of m. 

This hypothesis is natural in view of Conjecture 4.3, which says that the coefficients 
fiirn) in (7.13) are given for m > £ + 1 by polynomials Fi{m) = tim + (f)(, of degree 1 in 
m; so it is natural to expect that J-mio) will likewise satisfy a bound that is linear in m. 
The trouble is that the coefficients /^(m) ior £ > m are completely uncontrolled, so we 
cannot deduce Conjecture 7.3 from Conjecture 4.3. It is, nevertheless, a natural extension 
of Conjecture 4.3. 

Let us now prove that Conjecture 7.3 implies Conjecture 7.2: 

Proof of Conjecture 7.2, given Conjecture 7.3. Using the hypothesis (7.15) 
together with the Cauchy integral formula in the variable z — 1/q, we conclude that 

\T^{q)\ < whenever 1^1 > 2Qi . (7.16) 

We now use the fundamental formula (7.7) determining the curve B^, and we consider 
Im J^m{q)/m as a function of argg for each fixed |g|: 

Ue^Tgq;\q\) = Im . (7.17) 

m 

It follows from (7.16) that 

/^(argg;|g|) < 4C^ whenever |g| > 2(^1 (7-18) 



and hence that 

fLi^^gq; \q\)\ < \ whenever |g| >Q2 = max{8CQi,2Qi) . (7.19) 

Now (7.7) amounts to solving the equation 

/i^(argg) = y (7.20) 
for a set of values of y spaced by 7r/m, where 

hm{x) = X + fm{x; \q\) (7.21) 
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is a continuously differentiable (in fact infinitely differentiable) function satisfying h'^{x) > 
1/2 (and also /^^(a;) < 3/2, but we do not need this) whenever \q\ > Q2- It follows that 
the spacing between consecutive solutions of (7.20) cannot exceed 27r/m. Taking m — >■ oo, 
we conclude that the union of the curves Bm is dense in the region \q\ > Q2. It easily 
follows from this and the Beraha-Kahane- Weiss theorem that the union of the chromatic 
roots of Sm,n as m, n — >■ OO is dense in the region \q\ > Q2 well. □ 

7.3 Size of the discontinuity in the derivative of the free energy 

It is curious that the branches at large |g| of the curves Bm are caused by very small 
differences between the two eigenvalues — namely, |A*- — A*,^^! ^ 1 compared to |A*-| ^ 
|A*,^| ~ \q\'^ — which moreover become irrelevant in the limit m — >■ 00. 

In Figure 16 we plot the real part of the free energy 

fm{q) = max[/=(g),/^(g)] = -max[logA*-(g),logA*,^(g)] (7.22) 

for a semi- infinite strip of width m(l<m<6)asa function of Req for fixed Imq — 3 
or 4. The sohd dots show the points of discontinuity in the derivative of the free energy, 
which arise from the transition between the dominance of Xi,^= and A^^^. The discontinuity 
is small already for m = 1 and Img = 3, and it decreases as m and Img grow; indeed, it 
is essentially invisible for all m > 2. 

This should be contrasted with the behavior at the phase transition between the "dis- 
ordered" phase (at large \q\) and the "ordered phase" (at small \q\), which are separated 
by a curve lying close to BcxD(sq) = hm B^^ (see Section 6.7). At this transition — which 

occurs also for cylindrical and free boundary conditions, and which is the fundamental 
phase transition for this model — the discontinuities in the derivative of the free energy 
do not disappear as m ^ 00. This can be seen clearly in Figure 17, which is analogous 
to Figure 16 but for smaller values of Imq (namely, 0, 0.25, 0.5 and 1). 

We conclude that the outward branches (tending to 5 = 00) of the curves Bm do 
not correspond to phase transitions of the infinite- volume free energy. This is no surprise, 
because no such phase transitions are observed for free or cylindrical boundary conditions, 
and we believe that the infinite- volume free energy f{q) is the same for =, 7^, free and 
cylindrical boundary conditions (at least in the "disordered" phase, where we have checked 
it in the large-g' expansion through order cf. Section 5.2). 

7.4 Largest magnitude of a chromatic root 

By Proposition 7.1 we know that, for each m > 1, the chromatic roots of Sm,n become 
dense as n — > 00 on a curve Bm that contains 2m outward branches extending to infinity. 
It is natural to ask: At what rate does this curve get filled out as n — >■ 00? In particular, 
what is the largest magnitude of a chromatic root of Sm,n^ 

For the generalized theta graphs 0^*'^'), Brown, Hickman, Sokal and Wagner [12, The- 
orem 1.3] proved that the chromatic roots are bounded in modulus by [1 + o(l)]p/ logp, 
uniformly in s and p, where o(l) denotes a constant C{p) that tends to zero as p — > 00. 
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They furthermore showed [12, Corollary 5.2] that this bound is sharp when s = 2. How- 
ever, when s > 2 we expect that this bound is far from sharp; in particular, we expect 
that the root of largest modulus grows as p — >■ oo only like (p/ logpy^^^'^\ 

For our graphs S^^n, the only thing we know rigorously is that the chromatic roots 
are bounded in modulus by Cmax(n, 4), where 

and W is the Lambert W function [17], i.e. the inverse function to a; i— )> xe^; this follows 
from a result of Fernandez and Procacci [26] bounding the chromatic roots of arbitrary 
graphs in terms of the graph's maximum degree (i.e., the maximum number of neighbors 
of any vertex). However, by analogy with the generalized theta graphs it is natural to 
expect that the chromatic roots of the graphs Sm,n in fact obey an upper bound that is 
sublinear in n, e.g. n/logn or even (n/ logn)^^™- [since s corresponds to m + 1]. 

It is an interesting open problem to verify (or refute) these conjectures, first for the 
generalized theta graphs and then for the family Sm,n- 



8 Discussion 



In this section we make some final comments and some suggestions for future research. 



8.1 Comparison with the work of Jacobsen, Saleur and Dubail 

In a series of recent papers [20,39-41], Jacobsen, Saleur and Dubail have studied 
models of densely packed self-avoiding loops on the annulus in which loops touching one 
or both rims of the annulus receive different weights from "bulk" loops. By a generalization 
of the Baxter-Kelland-Wu [4] mapping, such models can be mapped onto Potts models 
in which spins lying on the rims of the annulus lie in specified subsets of {1, . . . ,q}, or 
equivalently Fortuin-Kasteleyn models in which clusters touching one or both rims of the 
annulus receive different weights from "bulk" clusters. 

These latter models can be defined in a very general way as follows [20]: Let G = (V, E) 
be a finite graph, and let Vi,V2 Q V be disjoint subsets of V. Fix subsets 81,82 Q 
{1, . . . ,q}, and let us define a Potts model in which spins at a site in Vi are required to 
lie in {i = 1,2): 

Z^'^''%q,8^,82,v) = J2 l[ [l + vM^i,aj)]. (8.1) 

a: V ^ {l,2,...,q} e=ij&E 
<r[Vi] C Si 
alV2] C S2 



^''Earlier, Sokal [62] had obtained the same result with the weaker constant 7.963907. The identification 
of the Fcrnandez-Procacci constant in terms of the Lambert W function is due to Jackson, Procacci and 
Sokal [31]. 
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Likewise, let us define a Fortuin-Kasteleyn model in which clusters touching Vi or V2 or 
both receive difi^erent weights from "bulk" clusters: 

Zaiq, gi, q., qu, v) = J] q'^^^^'^^^^'.^^^VA'^^^ J] , (8-2) 

ACE eeA 

where A;o(^), ^i(^), ^2(^), ^i2(^) denote, respectively, the number of connected compo- 
nents in the subgraph {V, A) that intersect neither Vi nor V2, Vi but not V2, V2 but not 
Vi, or both Vi and V2. It is then easy to see that 

ZP°"«(g,5i,52,v) = ZG{q,qi,q2,qi2,v) (8.3) 

provided that we identify 

gi = l^il, ?2 = 1^2!, gi2 = \SinS2\ . (8.4) 

Our partition functions Zs=^{q, v) and {q, v) are, up to trivial prefactors, special 
cases of this construction. Indeed, let G be the usual m x n square lattice with free 
boundary conditions (m > 2), and let Vi (resp. V2) consist of the sites in the leftmost 
(resp. rightmost) column. Then it is easy to see that 

^s=,„(g,v) = qZG{q,q-l,q-l,q-l,^r) (8.5a) 
^st„(^'^) = ?(?-l)^G(g,g-l,g-l,g-2,v) (8.5b) 

Dubail, Jacobsen and Saleur (DJS) have used the loop-model representation associated 
to the Fortuin-Kasteleyn model (8.2) to obtain the conformal properties (central charge, 
critical exponents, exact continuum-limit partition functions, etc.) both for the usual 
critical point n = [20] and for the Berker-Kadanoff phase n = — [41], where n is 
the bulk loop weight. This analysis is relevant for values of q lying inside the limiting 
curve Soo(sq). By contrast, in the present paper we have focused on the behavior at 
large \q\, i.e. on the disordered (and non-critical) phase lying outside the limiting curve 

There are, however, some points of contact between our work and that of DJS. In 
particular, our transfer-matrix eigenvalues should correspond to the sector with no non- 
contractible clusters in the work of DJS (who used periodic longitudinal boundary condi- 
tions) . 

8.2 Density of zeros 

Let us go back to the question raised at the beginning of this paper: When does the 

accumulation of partition-function zeros at some parameter value (or on some curve in 
the complex plane of some parameter) signal a nonanalyticity of the infinite-volume free 
energy? The answer is straightforward: when the zeros have a nonzero density per unit 
volume in the infinite- volume limit. Thus, in the trivial example Zn{x) = x for a system in 

^^The analysis of the Berker-Kadanoff phase in [41] is actually restricted to the one-boundary case 
92 = 5, (?12 = 91, but it could presumably be generalized. 
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volume n, the zero at a; = has multiphcity 1, so the density of zeros is 1/n, which tends 
to zero as n — )■ oo; hence the infinite- volume free energy does not exhibit nonanalyticity 
X — 0. Similarly for an isolated limiting point in the Beraha-Kahane-Weiss theorem: 
it has fixed multiplicity, independent of the system length n, so that the density of zeros 
again tends to zero as n — >■ oo. 

It would therefore be very interesting to try to understand analytically the density 
of zeros in the various graph families considered in this paper: starting with the bi-fans 
Pn + K2 and bipyramids C„ + K2 (Section 2.4), continuing with the generahzed theta 
graphs ©(*'P) (Section 2.5), and finishing with the graphs Sjn,n- For the first three families 
one might hope to find exact expressions for the density of zeros as n or p tends to infinity; 
furthermore, for the generalized theta graphs one might be able to study the subsequent 
limit s ^ 00. For Sm,n it seems unlikely that one could obtain exact expressions (except 
for very small m); but one might hope at least to obtain expansions in powers of q~^. 

8.3 Some concluding remarks 

In this paper we have studied the graphs Sm,n that are obtained from an m x n square 
grid (m columns, n rows, free boundary conditions in both directions) by adjoining one 
new vertex adjacent to all the sites in the leftmost column and a second new vertex 
adjacent to all the sites in the rightmost column. By similar transfer-matrix methods 
one can study the graphs Tm,n that are obtained analogously from a m x n grid of the 
triangular lattice. (The triangular lattice is most easily viewed as a square lattice with a 
SW-NE diagonal edge added to each face.) All of these graphs are 3-connected and planar, 
and we expect the square and triangular families to have a similar qualitative behavior, 
though the details of the limiting curves will differ [34,38,53]. The triangular lattice 
has the advantage that there exists a precise conjecture for at least the outer boundary 
of the limiting curve B^o for ordinary (e.g., free or cylindrical) boundary conditions: see 
Baxter [2,3] for this conjecture and the remarkable Bethe-Ansatz computation supporting 
it, and see [38, especially Section 6] for further discussion. 

More interestingly, one can also study the graphs S^l^ (resp. that are obtained 

by this construction starting from an m x n grid with periodic boundary conditions in the 
longitudinal (n) direction. These are 4-connected planar graphs (resp. 4-connected plane 
triangulations) , and for m = 1 they reduce to the bipyramids Cn + K2. Their chromatic 
polynomials can once again be analyzed by transfer-matrix methods [35,36]. The transfer 
matrices are larger than in the free-longitudinal-boundary-condition case, since they have 
to keep track of the connections among the sites of the top and bottom rows. On the other 
hand, the transfer matrix is block-triangular: firstly in terms of the number £ of disjoint 
paths connecting the top and bottom rows, and secondly in terms of the connectivity of 
the bottom row. Moreover, the amplitudes ak{q) are explicitly determinable polynomials 
in q. For this reason the computation is more involved than for free longitudinal boundary 
conditions, but not vastly computationally more demanding. See [35-37] for more details. 

Let us conclude with a vague idea. It may be possible to prove Conjecture 7.3 by 
using the polymer expansion [49] and exploiting the uniform (in m) convergence of the 
polymer expansion for sufficiently large j^l. 
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A Proof of Proposition 3.5 

In this appendix we prove Proposition 3.5 concerning the dominant diagonal entries 
in the transfer matrix. 

In Ref. [54, Section 3.2], we showed that the dominant diagonal entry tF(m) of the 
transfer matrix for a square-lattice strip of width m and free boundary conditions is 
equal to the partition function for a certain one-dimensional m-site polymer gas (with 
free boundary conditions). In this polymer model, each polymer of length £ > 1 gets a 
fugacity /i^ = v^~^{q + iv'), where v (resp. v') is the weight of all horizontal (resp. vertical) 
edges in the corresponding Potts model. We solved this problem by using the generating 
function ("grand partition function") 

oQ \]ir( \ 



m=l 

where 



^(^) = ^z^v^-\q + lv') 



l=\ 



+ 



\ — zv (1 — zvy 



(A.2) 



is the total weight for a single polymer of arbitrary size. It follows from (A.1)/(A.2) that 

We then expanded (A. 3) in powers of z to extract the coefficients tp{m). A similar 
analysis [54, Section 3.3] handled the case of cylindrical boundary conditions. 
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In this appendix we shall perform the analogous computations for the transfer matrices 
associated to the graph families S*^ „ and 5*^ „, culminating in a proof of Proposition 3.5. 

Let us start with the family As explained in Section 3.3 above, the boundary 

conditions for this family can be handled by considering m + 1 sites on a circle, where 
the site labelled is special. Consider first the action of H on the start vector eid. It 
generates 2™+^ terms, each of which corresponds to a partition V in which all the blocks 
are sequential sets of vertices on the (m + l)-cycle (we shall call these sets "polymers"). 
Each polymer of size I < m -\- 1 picks up a factor v^~^, while a polymer of size m + 1 
picks up a factor v'^^^ + (m + l)v'^ (the v'^^^ comes from the case in which all edges are 
occupied, while the (m + comes from the m + 1 cases in which all edges but one are 
occupied) . 

Consider next the action of V on a basis vector e-p corresponding to an arbitrary 
partition P = {Pi, . . . , P^}. If we are to end up with the partition eid, then for each block 
Pj we must either choose the detach operator Dj for all i e Pj (the last deletion gives a 
factor of q) or else choose the detach operator for all but one i G Pj and choose v'l for the 
last site (this can be done in \Pj\ ways). This means that the fugacity for any polymer of 
size 1 < £ < m not containing the special site 0, the fugacity is given by 

= v'-\q + ev'). (A.4) 

However, the fugacity for a polymer of length 1 < £ < m + 1 containing the special site 
is given by the formula^^ 

= forl<£<m 
^ \v"'{v + m+l) for£ = m + l ^ ' 

because the vertical operator V does not include Do (therefore, for any block Pj containing 
the site 0, we are forced to choose the detach operator Dj for all i e Pj other than the 
special site e Pj). 

As we did for cylindrical boundary conditions in [54, Section 3.3], we can compute 
t={m) by using a simple recursion relating our case to that of free boundary conditions: 

m 

t={m) = J2 k^t^ ^F(m + l-k) + ■ (A-6) 

fc=i 

To prove (A. 6), let A; > 1 the size of the polymer containing site 0. If k < m, there are k 
ways of placing this polymer such that the site belongs to it, with fugacity v^'' = v''~^ 
for each placement; and for the rest of ring, the total weight of all admissible polymer 
configurations is tF{m + l — k). Finally, if A; = m + 1, there is only one way of placing the 
polymer, and it receives a weight I'^^i- This proves (A. 6). 



^^Let us remark that the analogous definition [54, eq. (3.51)] suffers from an unfortunate notational 
ambiguity: /t^ should have been written as because the weight in question depends on m as well 

as £. The same superscript should appear also in [54, eqs. (3.52) and (3.53)]. This notational clarity 
is important because, although at this stage we are considering one fixed value of m, we will soon be 
summing over m. 
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In order to compute explicitly the t={m), we introduce the generating function 

oo 

$^(^) = ^;2"t=(m) . (A. 7) 



m=l 



The upper limit on the sum in (A. 6) can be replaced by oo, provided that wc define 
^f(^) = for £ < [which is anyway implicit in the definition (A. 3)]. Multiplying both 
sides of (A. 6) by z"^ and summing over m, we arrive at the equation 

(1 — vzy 1 — vz (1 — vzy 

When V = v' = —1, we obtain the final formula 

By expanding this function in powers of z, we have checked that it agrees with the 
dominant diagonal elements t={m) for m < 15. 

Let us next consider the family S^^^. Defining the generating function 

oo 



m=l 



it follows easily from (A. 7) and (3.13) that 



v'^z 



^A^) = $=(^) - — — . (A.ll) 

-L U /C 

When V — v' = —1, we obtain the final formula 

<^,(z) = -1 + — — . (A.12) 

' 1- {q-3)z- {q-l)z^ ^ ' 

By expanding this function in powers of z, we have checked that it agrees with the 
dominant diagonal elements t-^{m) for m < 15. 

It is now straightforward to extract from (A.7)/(A.9) and (A.10)/(A.12) the coeffi- 
cients t={m) and t^(m). In view of (3.13)/ (3.14), it suffices to study one of the two; and 
it turns out that t^{m) takes a shghtly simpler form [compare (3.24a) with (3.24b)]. Let 
us use the notation [z"^]P{z) to denote the coefficient of z^ in a polynomial or formal 
power series. Using the identity [54, eq. (3.21a)] 

1^-113^= E ('7 (A.13) 

we get for m > 1 

tAm) = [z-\^A^) = J] r . (g-3)-2^(g-l)7 (A.14) 

j=o \ J / 
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which is manifestly a polynomial in q of degree m, as claimed in (3.24). 

The next goal is to compute the coefficients ak{m) arising in (3.24), where m > 1 and 
< A; < ni. Expanding the binomials in (A. 14), we have 

We now want to substitute the m-dependent upper index in the sum over j by something 
independent of m, e.g. by k. 

There are two nontrivial cases: a) If A; < [m/2j, then the second binomial vanishes 
whenever k + i — 2j < 0, and the third binomial is non- vanishing only if j > i. Therefore 
for j > k and j > i we have that k + i — 2j<k + i — k — i = 0. So all these terms vanish, 
b) If A; > [m/2j, then the first binomial does not vanish only when < j < [m/2j or 
when j > m + 1. As we are adding terms with |_m/2j + l<j<k — l<'m — 1, none of 
them give rise to a non-vanishing contribution. Therefore, we can re-write (A. 15) as 

where the independent variable m docs not appear in the summation limits. After some 
straightforward but lengthy algebra we can rewrite the above formula in a more compact 
form: 

«^M^Err)M.'EK'"r^)G_:_j. 

It is clear from (A. 17) that afe(m) is for fixed A; > 0, the restriction to integers m > 
max(l. A;) of a polynomial in m of degree at most k, asm appears only in the upper index 
of the binomials and 

m\ m{m — l){m — 2) • ■ • {m — j -\- 1) fA 18) 

is a polynomial in m of degree j. [Here we use Knuth's [30] notation for falling powers: 
xi — x{x — 1) • • • {x — j -\- 1)] 

To see that the degree of ak{m) is exactly k, let us extract the term of order m'^: 

\m!']akim) = V \,7 = ^ ^ . (A.19) 

This equation also imphes that the coefficient akfi — 1 [cf. (3.30a)] for all A; > 0. 
This completes the proof of Proposition 3.5. 

Just to be safe, we have checked for < A; < m < 30 that (A. 17) agrees with the 
expansion of (A. 12). 
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B Dimension of the transfer matrices 



In this appendix we compute the dimensionahties N={m) and N^{m) of the transfer 
matrices for = and 7^ boundary conditions, proving the formulae stated in Theorems 3.6 
and 3.7. As a corollary we obtain the dimensionality dimT"(m) = SqFree(m + 2) = 
N^{m) + N=[m) of the full symmetry- reduced transfer matrix, confirming the formula 
(2.17) and proving the alternative formulae (3.40a,b,c). 

B.l Catalan, Motzkin and Riordan numbers 

The number of non-crossing partitions of {1, ... , m} is given by the Catalan number 



[66] 

r7_. = 

m\{m + l)\ m + l\ m 
The generating function for the Catalan numbers is 

C{z) = ^ C^z- = ^"^^^ . (B.2) 

m=0 

The number of non-crossing non-nearest-neighbor partitions of the set {1, 2, . . . , m} is 
given by the Motkzin number M^_i, where [9]^*^ 

Lm/2J , . 

Mm - Y: (^j^]c,. (B.3) 



fe=0 



If we consider the same vertex set on a circle, then the number of non-crossing non- 
nearest-neighbor partitions is given by the numbers [53] 

, f 1 for m = 1 /r, A\ 

^--{Rr. for m > 2 ^^-4) 

where the Rm are the Riordan numbers (or Motzkin alternating sums) defined by [9]^^ 

1 for m = 

= { m-l (B.5) 

E (-l)'""*^"^Mfe for m > 1 
fe=i 

In the following we will need the linear recursion for the Riordan numbers 

Rm — —Rm-l + Mm-1 + Sm,0 (B-6) 



Warning: Several references use the notation to„ to denote what we caU M„; and one reference [19] 
writes M„ to denote a different sequence. 

^^Equation (3.59) of [53] erroneously runs the sum down to A; = 0. 
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and the generating functions for M^, Rm and dm'- 



m=0 



^ ^ 2^(1 + ^) 



m=0 



(B.8) 
(B.9) 



, ^ N 1 - z + _ ^(1 + ^)(i -3^) 
= = R{z)-l + z = 241! . 

m=l ^ ' 

The numbers for 1 < m < 16 are displayed in Table 14. 

B.2 Partitions on a circle modulo reflection 

Our goal is to compute the number of equivalence classes modulo reflection of non- 
crossing non-nearest-neighbor (or "ncnnn" for short) partitions of the set {1,2,..., m) on 
a circle. However, on a circle there are two distinct types of reflection: 

• Reflection V.i with respect to an axis going halfway between a pair of neighboring 
sites (which we shall fix to be 1 and m) 

• Reflection T^-s with respect to an axis going through a site (which we shall fix to be 
the site 1) 

Therefore, we define N^irn) [resp. Nj,{ra)\ to be the number of equivalence classes modulo 
a refiection of type 1^2 (resp. 7^3) of ncnnn partitions of the set {1,2,..., m} on a circle. 
Of course, the two types of refiection coincide when m is odd, so that 

iV2(2£ + l) = iV3(2£+l) for£ = 0, 1,2,... . (B.IO) 

For m even we have to solve the two problems separately. 
The first step is to define 

X{m) — #ncnnn partitions of {1, . . . , m} on a circle 

that are invariant under a refiection of type 1^2 ■ (B-H) 

The values of X{m) for 1 < m < 16 are displayed in Table 14. We furthermore denote by 
X(m) the set of all 7?.2-invariant ncnnn partitions of {1, ... , m} on a circle. The first sets 
X{m) are easy to write down. Using the delta-function shorthand (2.13) for partitions, we 
find that X{m) = {1} for m = 1, 2, 3, 4, X(5) = {1, ^24}, ^(6) = {1, ^25, Srdm}: X{7) = 

{1, ^26, ^35, ^26^35, ^13^57, ^246}, and X(8) = {1, 627, S36, 527^36, 513^68, ^14^58, 524557}- 

Because each equivalence class modulo refiection contains either one or two elements, 
according as those elements are or are not invariant under refiection, we have 

N2{m) = ^[X{m) + dm] . (B.12) 
It turns out that the numbers X(m) give us not only A^2 but also N^: 
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Lemma B.l The number N^^rn) of equivalence classes modulo a reflection of type TZ^ of 
ncnnn partitions of the set {1,2, ... , m} on a circle is given by 



Proof. The equality for odd m is trivial from (B.10)/(B.12). For even m wc start by 
defining the quantity 



which is equal to the number of ncnnn partitions that are invariant under a refiection of 
type 7^3. Let us now show that Y{2i) = X{2i + 2) for a\\ £ > 1. In this proof we shall 
denote by Y{m) the set of all 7^3-invariant ncnnn partitions of {1, . . . ,m} on a circle. 

The idea is to show that there is a bijection between the set Y{2£) of T^s-invariant 
ncnnn partitions of {1, 2, ...,£,£ + 1, ... , 2£} on a circle and the set X{2£ + 2) of TZ2- 
invariant ncnnn partitions of {1,2, ...,£,£ + 1, (£ + 1)', £ + 2, . . . ,2£, 1'} on a circle. The 
trick is to split the vertices 1 and £ + 1 into two pairs of nearest-neighbor vertices 1, 1' 
and £ + 1,{£ + 1)', respectively. 

This bijection is clear between the subset of Y{2£) with 1 and £ + 1 singletons, and 
the subset of X{2£ + 2) for which 1, 1', {£+l),{£+ 1)' are all singletons. Using the delta- 
function shorthand (2.13), those partitions belonging to the above subsets have exactly 
the same expression (e.g., for £ — 3, the ncnnn partitions ^2,6 and 83^5 belong to these 
subsets). 

There is also a bijection between the subset of Y{2£) with £ + 1 a. singleton and 1 
joined to some T^a-symmetric block B U B' with B C {3,A, . . . ,£} and B' — TZ^ B, and 
the subset of X{2£ -\- 2) with both {£ -\- 1) and {£ -\- 1)' singletons, and the sites 1 and 1' 
joined to the blocks B and B' , respectively. This property is due to the facts that 1 and 1' 
[resp. £+1 and {£+ 1)'] are nearest-neighbor vertices, and the partitions are non-crossing. 
For instance, if ^ = 3, the partition 5i 3 5 in 1^(6) corresponds uniquely to the partition 
^1,3^1', 5 in -^(8)) and vice versa. A similar bijection can be shown when 1 is a singleton, 
but ^ -|- 1 is joined to some T^-a-symmetric block. 

Finally, there is a bijection between the subset of Y{2£) with both 1 and £ + 1 joined 
to some 7^3-symmetric blocks Bi U B[ and B2 U B2 (that might be the same), and the 
subset of X{2£ + 2) characterized by the blocks {1} U Bi, {!'} U B[, {£ + 1] D B2, and 
{{£ + 1)'} U B'2. (If 1 and £ + 1 belong to the same block, the modifications are obvious.) 

As the above subsets are non-intersecting and contain all partitions of both Y{2£) and 
X{2£ + 2), we conclude that such bijection exists, and therefore, that Y{2i) — X{2i + 2), 
as claimed. □ 

Henceforth we consider only refiections of type TZ2- Our goal is to obtain a closed for- 
mula for the quantity X{m). Our results can be summarized in the following proposition: 

Proposition B.2 The quantity X{m) is given by 




N2{m) = \ [X{m) + d^] for m odd 
^ \Xim -I- 2) -I- drr,] for m even 



(B.13) 



Y{m) = 2N3{m)-d, 



(B.14) 




if m — 1 



X{m) 



) 



ifm > 2 



(B.15) 
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The generating function for this sequence is 



oo 



1 -z'^ + l 



G{z) = $^X(m) 



z' 



,m 



2Z ^(l + ^2)(l_3^2) 



+ 2^2-1 



(B.16) 



m=l 



In what follows, it is convenient to change the notation for the vertices: When m = 2i 
is even, we renumber the vertices 1, . . . , £, £ + 1, . . . , 2£ as 1, ...,£,£',... ,1', so that the 
vertices j and j' (1 < j < i) are transformed into each other under a reflection 7^2- 
When m = 2£ + 1 is odd, we renumber the vertices 1, . . . , £, £ + 1, £ + 2, . . . , 2£ + 1 as 
1, ...,£,£+ 1,£', ... ,1', so that the vertices j and j' (1 < J < £) are again related by 
reflection, while the vertex £ + 1 is invariant under reflection. See Figure 18 for a few 
examples. 

Proof of Proposition B.2. The idea is to obtain a recursion for the X{m), and 
then solve this recursion. First, we want to obtain X{2m + 1) from X(2m). We use the 
notation i -H" j to denote that i and j belong to the same block of the partition being 
discussed. We observe that for a partition in X{2m), we have 1 <H- m if and only if 
1' ■<->■ m' (since the partition is 7?.2-invariant) . We then split the set X{2m) into three 
subsets, as follows: 

a) Xa{2m) is the set of partitions in X{2m) such that there is no block containing 
both unprimed and primed vertices and such that 1 m; 

b) Xfj{2m) is the set of partitions in X{2m) such that there is no block containing 
both unprimed and primed vertices and such that 1 ^ m; 

c) Xc{2m) is the set of partitions in X{2m) such that there is at least one block 
containing both unprimed and primed vertices. [Note that any such block B is 
necessarily 7^2-invariant, for otherwise B would cross 7^2-B-] 

We denote by Xa{2m), Xh(2m) and Xc{2m) the number of elements within the correspond- 
ing subset. For instance, for X{8) we have the decomposition Xa{S) = {1, 61361131, 62a62'4'}, 
Xb{8) = {6i4,6r4'} and Xc{8) = {^22', fe, 522'533'}- 

Let us suppose that m > 3. Going from X{2m) to X{2m + 1) means that we have to 
add a new vertex labeled m + 1 between m and m'. Indeed, if e X(2m), then it also 
belongs to X{2m + 1) [with the new vertex m + 1 becoming a singleton]. Furthermore, 
to any partition V € Xa{2m) U Xc{2m) we can adjoin a factor 6mm' [i-e. contract the 
blocks containing m and m' into a single block] and get a partition in X{2m + 1) [with 
the new vertex m + 1 again becoming a singleton] ; please note that this is not possible if 
V e Xii{2m) since the new partition would include a nearest- neighbor connection 1 ^ 1'. 
Finally, if V E Xc{2m), we can generate a new partition in X{2m + 1) by joining the 
new vertex m -|- 1 to the closest block containing both primed and unprimed vertices. It 
is not hard to see that in this way we generate all the elements of X{2m + 1), and that 
we generate each one exactly once; therefore. 



X{2m + 1) 



X{2m) + Xa{2m) + 2X^{2m) . 



(B.17) 
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Now, Xa{2m) +Xh{2m) is just the total number of partitions in X{2m) with no block 
containing both unprimed and primed vertices. But this number is obviously equal to the 
number of nennn partitions of the set {1, 2, ... , m} with free boundary conditions, i.e. the 
Motzkin number M^_i: 

Xa{2m) + Xb(2m) = M^-i . (B.18) 

On the other hand, Xi,{2m) can be interpreted as the number of ncnnn partitions of 

{l,2,...,m} with free boundary conditions, such that m and 1 always belong to the 
same block. But this equivalent to the number of ncnnn partitions of{l,2,...,m — 1} 
on a circle, so we have 

Xb{2rn) = dm-i = Rm-i (B.19) 
since m > 3. Putting everything together we obtain the recursion 

X{2m + 1) = 3X{2m) - {Mm-i + Rm-i) for m > 3 (B.20) 

together with the initial conditions ^(1) = 1, -^^(3) = 1, X{b) — 2. It is interesting to 
note that the same recursion holds also for m = 1, 2, as can be checked from Tables 14 
and 1. Thus, our final result is 

X{2m + 1) = 3X(2m) - {M^-i + Rm-i) for m > 1 (B.21) 

with the initial condition X{1) = 1. 

The second half of the recursion can be obtained by considering how to obtain the 
partitions in X{2m + 2) by starting from those in X{2m + 1). In this case, the vertex 
m + 1 is split into a pair of nearest- neighbors m + 1 and (m + 1)'. Let us assume that 
m > 2. 

Consider a partition V G X{2m + 1) in which the block containing the vertex m + 1 is 
{m + 1} U 5 U B', where B contains only unprimed vertices and B' = contains only 
primed vertices (note that B could be the empty set). We can then create a partition 
V e X(2m + 2) by creating blocks {m + l}U B and {(m + 1)'} U B' and leaving all 
other blocks of V as is. The partition V has the property that, if one contracts the 
blocks containing m + 1 and (m + 1)' into a single block, one obtains V; moreover, V 
is the unique partition in X{2m + 2) with this property (it is not hard to sec that all 
other 7^2-iiivariant partitions with this property would contain either a nearest-neighbor 
connection or a crossing). Moreover, every partition in X(2m-|-2) is obtained in this way 
except those in which 1 m + 1 [and hence also 1' ^ {m + 1)']; these latter partitions 
cannot be obtained because contracting m-|- 1 to (m-|- 1)' would create a nearest-neighbor 
connection 1 ^ 1'. Now, the number of these extra partitions is equal to the number 
of ncnnn partitions of {1,2,..., m} on a circle, which is dm] and since m > 2, this also 
equals Rm [cf. (B.4)]. We therefore conclude that 

X{2m + 2) = X{2m + 1) + Rm for m > 2 (B.22) 

together with the initial conditions X(2) — X[4) = 1. It is interesting to note that the 
same recursion holds also for m = 1, as can be checked from Table 14. Thus, our final 
result is 

X{2m + 2) = X{2m + 1) i?^ for m > 1 (B.23) 
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with the initial condition X{2) = 1. 

By combining (B.21)/(B.23) and (B.6), we can obtain the following recurrences for 

x{my. 

X(2m) = 3X(2m - 2) - 2R^_2 for m > 2 (B.24a) 

X{2m + 1) = 3X{2m -1)-Rm + Rm-i for m > 2 (B.24b) 

with the initial conditions X{m) = 1 for m = 1, 2, 3. 

A solution of (B.24a) can be obtained by defining X{m) — X{2m-\-2) for m > 0. The 
recurrence (B.24a) can then be written in terms of the X as 

X(m) = 3X(m - 1) - 2Rm-i for m > 1 (B.25) 

with the initial condition X(Q) — 1. From this equation we obtain the corresponding 
generating function 

G{z) = yX{m)z- = iz^fM ^ ^ , (B.26) 

where R{z) is the generating function (B.8) for the Riordan numbers. 

Let us now observe that (B.26) coincides with the generating function of the central 
trinomial coefficients [16, p. 163] [61, sequence A002426], so we have 

X{m) = [f^jil + 1 + 1^)"" form>0. (B.27) 

By a double application of the binomial theorem we arrive at the expressions 

Note that many terms on the right-hand side of (B.28) vanish, as (™^'^) ~ whenever 
k > m/2. In order to lighten the notation, we will use the generic summation over 
nonnegative k, and only at the end we will make the upper bound explicit. 

The formula for X{2m) and its generating function G2{z) can be obtained from 
(B.28)/(B.26): 

fe>0 \ / \ / 

2 



G,{z) = Y.X{2m)z'- = z'G{z') = . (B.29b) 

^1 + 



A solution of (B.24b) can be obtained in a similar way by defining X{m) = X(2m + 1) 
for m > 1 (we will treat -'^(l) = 1 separately). The recursion (B.24b) can then be written 
as 

X{m) = 3X{m -1)-Rm + Rm-i for m > 2 (B.30) 
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with the initial condition X(l) = 1. From this equation we obtain the corresponding 
generating function 



k=l ^ 



(B.31) 



(1 + z)(l - 3z) + (1 - z)y/{l + z){l-3z) 

We now observe that (B.31) coincides with the generating function of the next-to- 
central trinomial coefficients [61, sequence A005717], so we have 

X(m) = + t + t^)"" = + t + t^r form>l. (B.32) 

By a double application of the binomial theorem we arrive at the expressions 

^w-i:(:)(:;f)-E(:)(r9 ^-^^ 

fc=0 ^ ^ ^ ' ^ k=0 ^ / ^ / 

The formula for X{2m + 1) and its generating function Gi{z) can be obtained from 
(B.33)/(B.31): 



X{2m + 1) = <^ ^ fm\ fm-k\ ^ ^ (B.34a) 




G^{z) = J]X(2m + l)^'"^+^ = z + 

m=0 

= z + = (B.34b) 

(1 + _ 3^2) + (1 _ z^)^{l + ^2)(l _ 3^2) 

Prom (B.29a)/(B.34a) it is easy to obtain (B.15). Prom (B.29b)/(B.34b) it is easy to 
obtain the formula (B.16) for the generating function G{z) — Gi{z) + G2{z). □ 

By Proposition B.2 together with (B.12)/(B.13), it is easy to compute closed expres- 
sions for N2{m) and A^3(m) as well as for their generating functions. We have 

{1 for m = 1 

f^V * J\k + Ilrntsodd]J 



and 



for m = 1 

Nsi^) = { ^^^fmf LfJ-^ ^ form>2 ^^-36) 

k + I[m is odd] 
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The corresponding generating functions are 

1 

N2{m)z'" = 

m=l 



F^iz) = J2N2{m)z^ = -[G{z) + D{z) 



1 



2z-' -z^ + 1 



^y{l + Z^){l-3z^) V 1 + Z 



l-3z 



+ 2z{2z - 1) 



(B.37) 



z — - 



m=l 



G,{z) + ^-l + D(z) 



1 /l-3z 
= + 



1 + 2^-^2 



(B.38) 



Note that 
so that 



4z\j 1 + z 4z^y{l + z^){l-3z^) 

Fsiz) - F2{z) = z'G{z') , (B.39) 
N3{2£ + 2) - N2{2£ + 2) = [t^-^]{l + t + t'^f for £ > . (B.40) 



B.3 The dimensions N^(m) and N=(m) 

We can now prove the theorems stated in Section 3.3. 

Proof of Theorems 3.6 and 3.7. As explained in Section 3.3, we have 

N^{m) = A^2(m + 2) 



(B.41a) 
(B.41b) 



The formulae (3.34)/(3.37) for N={m) and N^{m) then follow trivially from the formulae 
(B.35)/(B.36) for N2{m) and N^{m). The results for the generating functions follow 
likewise from the identities 



G={z) = -[F,iz)-z] 
z 



z — z 



(B.42a) 
(B.42b) 



□ 



Remark. We see from (B.IO) and (B.41) that 

N^{m) = N={m + 1) for m odd 

(cf. Table 1). 



(B.43) 
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Proof of Corollary 3.8. From Theorems 3.6 and 3.7 we have 
dimT"(m) = N=(m) + N^(m) 



m' — k 



1 / ^ ^ X 1 fm'\ f fm' - k\ 



1,. m'\ m' -k + l\ 



where we have used the relation (B.6) and the elementary recursion relation for the 
binomials; this proves (3.40a). 

('yY)j \ ( 777, /u \ 
ki\ k / summand on 

both sides; using the elementary recursion for binomials, we get 

Efm\fm-k + l\ fm\ fm-k + l\ ^-n, ,r\ 

Applying this with m replaced by m' shows the the equivalence of (3.40b) with (3.40a). 

Finally, by adding (3.40a) and (3.40b), dividing by 2, and using once again the recur- 
sion relation for binomials, we prove (3.40c). □ 



Proof of (2.17). We replace m + 2 by m in (3.40b), which entails replacing m' by 
m' — 1. Expanding the binomials into factorials yields (2.17). □ 



C Chromatic polynomials for the family Sra,n 

We denote by Sjn^n the graph obtained from the square- lattice strip of width m, length 
n and free boundary conditions in both directions by attaching one extra site to all the 
sites on the top row and another extra site to all the sites on the bottom row (see Figure 8). 
This graph is obviously isomorphic to the graph Sn,m\ we have merely rotated the picture 
by 90 degrees. In particular, the chromatic polynomials of Sm,n and Sn^m must coincide, 
i.e. 

Ps^M = Ps^,r.{Q) , (C.l) 

and we can use this as a check on the correctness of our computations. 

Here we shall develop the transfer- matrix formalism for the family Sm,m where m is 
fixed (and small) and n is arbitrary; this is complementary to the formalism developed in 
Section 3, which treats Sn-m for n fixed (and small) and m arbitrary. 

The graph is is simply the path P„+2, so its chromatic polynomial is 



76 



The chromatic polynomial of the graph S'2,n can be computed by noting that this 
graph is obtained by gluing together n — 1 squares C4 and two triangles C3 along edges. 
Its chromatic polynomial is therefore 

To handle the graphs Sm,n with m > 3 we need a different technique. The idea is 
to use the standard transfer matrices T"(mF) for square-lattice strips with free boundary 
conditions'^ but with different initial and final vectors to take account of the extra sites 
attached to the top and bottom rows. More specifically, we have 

= ^lr'{m,r-% , (C.4) 

where the initial vector f_|_ is given by 

f+ = HVecontr (C.5) 

where econtr is the basis vector for the partition with all sites in the same block, while the 
final vector a;+ is defined by [compare to (2.11)] 

u;^ep = 1)1^1. (C.6) 

The reason for this last formula is clear: the factor q comes from the extra vertex on top, 
and we obtain a factor oi q — 1 for each block of the partition V (instead of a factor of q) 
because each site on the top row is attached to the extra site, so that there are only q — 1 
allowed colors for those sites (rather than q). 

For instance, for m = 3 the chromatic polynomial of the graph S^^n is given in the 
basis B = {l,5i3} by'' 



with 

g' - + lOg - 8 - 4g + 5 



T"(3.) = (^^ - ^ \-2 )■ 

For m — A, the corresponding chromatic polynomial is given in the basis B — {l,5i3 + 
Su} by 

{q-2f ^ ^ 

P^^^ = g(g-l)(g-2) I 2(g-2) I .T"(4pr^. I 1 | (C.9) 

q-3 




^^See Section 2.2 above, and see [53] for more details. Since v = —1, the horizontal transfer matrix 
H is a projection, so we can use the modified transfer matrix T' = HVH and the modified basis vectors 
f-p = He-p. Moreover, since the graphs Sm,n are invariant under reflection in the central vertical axis, we 
can use the symmetry-reduced transfer matrix T". These transfer matrices are precisely the ones that 
can be found in [53, Section 5] (but denoted there T rather than T"). 

^^Please recall from Section 6 that we are using a shorthand notation for basis elements: we write 
the basis as a collection of vectors e-p, but it should be understood that the basis vectors are actually 
fp = Hep. This convention applies throughout this Appendix. 
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with 



jq^ -Iq^ + 2lq^ -2,2q + 21 2{q^ - Qq"^ + lAq ^ 12) q^ - Iq^ + 1% - 2Q\ 
T"(4f) = q-2 q^-Aq + b 3-g 

V -1 -2(g-2) q^-bq + 7 ) 

(C.IO) 

We have performed the corresponding computation also for widths m — 5 and 6, and we 
have checked that the identity (C.l) holds for 1 < m < 6 and 1 < n < 8. 

We can alternatively obtain the chromatic polynomial of the strip graph Sm,n by using 
a different technique, which is useful when dealing with more complicated endgraphs. Our 
starting graph consists of a strip together with attached endgraphs at top and bottom; 
each of these endgraphs consists of one or more extra vertices and several extra edges. 
The idea is to use the deletion-contraction identity on the extra edges of the endgraphs to 
get rid of the extra vertices. At the end, the chromatic polynomial Pg can be written as 
a linear combination of chromatic polynomials Ps^ „ associated to different strip graphs 
Sm,n,j, each of which is obtained from the mp x np grid by adding only extra edges to 
the top and/or bottom rows. (This approach is closely related to a method developed 
by Rocek, Shrock and Tsai [52].) These extra edges can be handled by inserting the 
appropriate operators Qij = I + VijJij at left and/or right in the standard expression 
(2.10) for the chromatic polynomial in terms of transfer matrices. There is, however, 
one complication: for m > 4, these extra edges can give rise to crossing partitions (for 
instance, if in a grid of width m = 4 we add an extra edge 13 on the bottom row and 
an extra edge 24 on the top row).^^ We must therefore work in a basis that includes also 
(non-nearest-neighbor) crossing partitions. 

For m = 3 we have 

Ps,Jq) = (g-2)^P3pxnp(?)-2(g-2)P53,„,.(g) + P53.„,(g) (C.ii) 

where 3f x np is a square-lattice strip of width m = 3, length n and free boundary 
conditions; S3^n,a is the chromatic polynomial of the graph obtained by adding to 3f x n-p 
an extra edge joining the sites 1 and 3 on the bottom row; and Ss^n,b is the graph obtained 
by adding to S3^n,a one further edge joining sites 1 and 3 on the top row. The chromatic 
polynomials of these three graphs can then be obtained by a transfer-matrix formalism. 
For the graph 3f x np, we can use the standard method discussed in Section 2.2. For 
the other two graphs, we obtain the corresponding chromatic polynomial by inserting the 
operators Q13 appropriately. Thus, the chromatic polynomial for Sa^n can be written as 

Pssjl) = (g-2)2p3^,„^(g)-2(g-2)a;T-T"(3F)"-^-Qi3-fid 

+u;T • Qi3 • T"(3f)"-' • Qi3 • fid (C.12) 

^"^Notice that this can happen even though the graph with the extra edges is planar (as it is in the 
just-mentioned example). To avoid crossing partitions we need to be able to draw the graph in a planar 
fashion within the m x n box. 
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where T"(3f) is given by (C.8). After some algebra we find 

T 



53,71 



q{q-l){q-2){{q-2) 



q-1 



T"(3p) 



n— 1 



+ 



3-2q 
-2 



T 



T"(3p) 



n-l 



(CIS) 



For m = 4 the situation is similar but more involved algebraically. We shall work 
in the basis B = {l,5i3 + ^24,614,613624}, where the last basis element corresponds to 
the crossing partition. After application of the deletion-contraction identity, and some 
algebra, we find 

Ps^Jq) = (g-2)2p4pxnp(g)-2(g-l)u.TT"(4Fr-iQi4fid 

-2iq - l)^a;TT"(4F)"-i(Qi3 + Q24)fid + u;^Qi4T"(4F)"-^Qi4fid 
+2a;TQi4T"(4F)"-^(Qi3 + Q24)fid 



with the transfer matrix 



T"(4f) 



24)T"(4f)"- 




13 + Q24)fid 


( 




q^-hq + 8\ 


T"(4f) 










-1 


^0 





1 / 



(C.14) 



(C.15) 



where the upper- left 3x3 submatrix is given by (C.IO). Notice that we have symmetrized 
the operators Q13 and Q24 in order to use the symmetrized connectivity basis B as above 
(i.e., classes of non-nearest-neighbor connectivities that are invariant under reflections 
with respect to the center of the strip). The definition (2.11) of the left vector uj continues 
to be valid even in the presence of crossing partitions. After some algebra we find 



M) 



q{q-l){{q-r 



f q^-3q + 3\ 
2{q-l) 

1 



2{q-l) 
q-2 
1 

T 



+ 2(g-2) 



T"(4f)"-' 
/ 5 \ 



fq-7\ 

2 
2 

V J 



T"(4f) 



n-l 



1 
1 

V / 



T"(4f) 



n-l 



-2 

-1 

V / 

/ 2 \ 

-1 


V / 



(C.16) 



Note that in this final expression the crossing partitions play no role, because all three 
right vectors have zero in the last entry, and the transfer matrix (C.15) has zeros in 
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the last row in all columns other than the last. So we could have restricted attention 
to non-crossing partitions, after all! We suspect that a cancellation of this type occurs 
whenever the graph is planar, and that a transfer-matrix description using only non- 
crossing partitions can be obtained by adding extra layers at top and bottom together 
with suitable v = oo vertical and/or horizontal edges. 

With this method we have checked the identity (C.l) for 1 < m < 4 and 1 < n < 8. 
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Table 1: Dimension of the transfer matrix. For each strip width m we give the number 

Bm of all partitions, the number Cm of non-crossing partitions, the number M^-i of non- 
crossing non-nearest-neighbor partitions with free boundary conditions, and the number 
SqFree(m) of equivalence classes of non-crossing non-nearest-neighbor partitions modulo 
reflection in the center of the strip. The dimension of the transfer matrix T"(m) associated 
to the strip Sm,n is SqPree(m -|- 2). We also show the dimensionalities N^{m) and N={m) 
of the transfer matrices T'A(m) and T^(m), respectively. 
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Table 2: Coefficients fc^ (m) of the large-g expansion of the dominant eigenvalue A*_^. For each 1 < m < 11, we include all 
coefficients &f (m) up to A; = 15. For the whole data set up to A; = 40, see the Mathematica file data_Neq.m included in the 
on-line version of the paper at arXiv.org. Those data points below the staircase-like line satisfy m > [cf. (4.16)]. 
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Table 3: Coefficients kcf^m) of the large-g expansion of log(g ™Af ) where A^^^ is the dominant eigenvalue for 7^ boundary 
conditions. For each 1 < m < 11, we include all coefficients c^{m) up to /c = 15. For the whole data set up to /c = 40, see 
the Mathematica file data_Neq.m included in the on-hne version of the paper at arXiv.org. Those data points below the 
staircase-like line satisfy m > mt^^^lk) [cf. (4.16)]. 
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Table 4: Polynomials Cfc(m), Ek{m) and Fk{m) representing for m > + 1 the coeffi- 
cients c}^,{m) [cf. (4.3)], efc(m) [cf. (4.13)] and fk{^) [cf. (4.14)], respectively. All these 
polynomials are of degree 1 in m, and their coefficients multiplied by k are always integers. 
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Table 5: Coefficients A^(m) for 7^ boundary conditions [cf. (4.21)] for 1 < A; < 15 and 1 < m < 11. Those data points below the 
lower staircase-like line satisfy m > rn^^^^^k) [cf. (4.16)] and are therefore zero. Those data points between the two staircase-like 
lines can be fitted to a polynomial Ansatz [cf. (4.23)-(4.27)]. 
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Table 6: Coefficients b^{m) of the large-g expansion of the dominant eigenvalue A*-. For each 1 < m < 12, we include all 
coefficients b^{m) up to A; = 15. For the whole data set up to A; = 40, see the Mathematica ffie data_Eq.in included in the 
on-line version of the paper at arXiv.org. Those data points below the staircase-hke hne satisfy m > [cf. (4.28)]. 
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Table 7: Coefficients kc^(m) of the large-g expansion of log(g^"' A^) where A^,= is the dominant eigenvalue for = boundary 
conditions. For each 1 < m < 12, we include all coefficients c^irri) up to k = 15. For the whole data set up to A; = 40, see the 
Mathematica file data_Eq.m included in the on-line version of the paper at arXiv.org. Those data points below the staircase-like 
line satisfy m > m^i^{k) [cf. (4.28)]. 
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Table 8: Coefficients Ar(m) for = boundary conditions [cf. (4.21)] for 1 < A; < 15 and 1 < m < 12. Those data points below the 
lower staircase-like line satisfy m > m^i^{k) [cf. (4.28)] and are therefore zero. Those data points between the two staircase-hke 
lines can be fitted to a polynomial Ansatz [cf. (4.29d)-(4.31)]. 
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Table 9: Coefficients di{m) [cf. (4.11)] for 1 < £ < 13 and 1 < m < 11. For the whole data set up to £ = 20, see the Mathematica 
file data_Dif f .m included in the on-line version of the paper at arXiv.org. Those data points below the lower staircase-like line 
satisfy m > i + I. 
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Table 10: Coefficients iei{m) [cf. (4.13)] for 1 < £ < 13 and 1 < m < 11. For the whole data set up to £ = 20, see the 
Mathematica file data_Dif f .m included in the on-line version of the paper at arXiv.org. Those data points below the lower 
stair case- like line satisfy m> £ + 1. 
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Table 11: Coefficients Af\m) [cf. (4.34)] for 1 < £ < 13 and 1 < m < 11. For the whole data set up to £ = 20, see the 
Mathematica file data_Diff .m included in the on-line version of the paper at arXiv.org. Those data points below the lower 
staircase- hke line satisfy m > £ + 1 and are therefore zero. Those data points between the two staircase-hke fines can be fitted 
to a polynomial Ansatz [cf. (4.35)]. 
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Table 12: Coefficients /i;/^(m) [cf. (4.14)] for 1 < £ < 12 and 1 < m < 11. For the whole data set up to i = 20, see the 
Mathematica file data_Diff .m included in the on-line version of the paper at arXiv.org. Those data points below the lower 
staircase- hke hne satisfy m> i + 
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Table 13: Coefficients Ay^(m) [cf. (4.37)] for 1 < £ < 12 and 1 < m < 11. For the whole data set up to £ = 20, see the 
Mathematica ffie data_Diff .m included in the on-line version of the paper at arXiv.org. Those data points below the lower 
staircasc-like line satisfy m > f + 1 and are therefore zero. Those data points between the two staircase-like lines can be fitted 
to a polynomial Ansatz [cf. (4.38)]. 
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Table 14: Counts of certain classes of non-crossing non-nearest-neighbor (ncnnn) parti- 
tions of the set {1, 2, ... , m} on a circle. For each m, we show the total number d„i of 
ncnnn partitions (which equals the Riordan number i?^ when m > 2); the number X{m) 
of ncnnn partitions that are invariant under reflection with respect an axis going between 
vertices 1 and m; the number N2{m) of equivalence classes of ncnnn partitions modulo 
reflection with respect an axis going between vertices 1 and m; and the number N^lm) 
of equivalence classes of ncnnn partitions modulo reflection with respect an axis going 
through vertex 1. Note that -Ro = 1 and Ri — 0. 
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12 3 4 



(a) (b) 

Figure 1: (a) Square- lattice strip 5*3,6 of width m = 3 and length n = 6 with two extra 
sites, (b) Alternate representation of S^^ in terms of a square-lattice strip of width 
m + 2 = 5, length n = 6 and free boundary conditions, with edge weights v — > +oo on 
the dashed edges. 
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Figure 2: Generalized theta graph B^^'^) formed by p = 5 chains in parallel, each consisting 
in s = 3 edges in series. 
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Figure 4: Chromatic roots for the bi-fan Pioo + K2 = Si, 100 (red •). We also show the 
limiting curve Bi, and the isolated limiting points (x). 
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Chromatic zeros of Limiting curves 0^^'^^: s=5,10,15 




Figure 5: Chromatic roots for the generahzed theta graphs Q^^'P\ (a) Chromatic roots 
for the graphs 0(^'P) with p = 25 (black n), p = 100 (red •) and p = 400 (blue A), and 
the limiting curve C5 [cf. (2.35)]. We also show the isolated limiting points for this family 
(black x). (b) Limiting curves Cs for s = 5 (black), s = 10 (red) and s = 15 (blue). The 
circle (depicted in green) is |g — 1| = 1. 
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(a) (b) 

Figure 6: (a) Square-lattice strip S^^ of width m = 4 and length n = 3 with the two extra 
sites (0 and 5) contracted. This contraction is indicated by a dashed line, (b) Alternate 
representation of S^^ in terms of a square-lattice strip of width m + 1 = 5, length n = 3 
and cylindrical boundary conditions, with edge weights v — t- +oo on the dashed edges. 
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(a) (b) 

Figure 7: (a) Square-lattice strip Sf^ of width m = 4 and length n = 3 with the two extra 
sites (0 and 5) joined by a f = — 1 edge, (b) Alternate representation of 5*1^3 in terms of a 
square-lattice strip of width m + 2 = 6, length n = 3 and cylindrical boundary conditions, 
with edge weights v — )■ +oo on the dashed edges. 
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Figure 8: Square-lattice strip 6*5^4 of width m = 5 and length n = A with extra sites at 
top and bottom. 



104 



Zeros sq lattice m = lp++ 



Zeros sq lattice m = 
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Figure 9: Limiting curves for square-lattice strips of width m = 1 with two extra sites. 
We also show the zeros for the strips S'i_io (black □) and S'1,20 (red o). We depict the 
isolated limiting points with the symbol x. The zeroth-, first-, second- and third-order 
asymptotic expansions for the outward branches [cf. (7.14)] are depicted as dashed orange, 
dotted red, dot-dashed magenta, and dot-dot-dot-dashed blue curves, respectively. Please 
note that for m = 1 the first- and second-order asymptotics coincide, as the term of order 
has the factor sin[(2A; — l)7r] = 0. 
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Zeros sq lattice m = 2^,+^ ^eros sq lattice m = 2- 




Figure 10: Limiting curves for square-lattice strips of width m = 2 with two extra sites. 
We also show the zeros for the strips 5*2,20 (black □) and 5*2,40 (red o). We depict the 
isolated limiting points with the symbol x. The zeroth-, first-, second- and third-order 
asymptotic expansions for the outward branches [cf. (7.14)] are depicted as dashed orange, 
dotted red, dot-dashed magenta, and dot-dot-dot-dashed blue curves, respectively. 
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Zeros sq lattice m = 3 
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Figure 11: Limiting curves for square-lattice strips of width m = 3 with two extra sites. 
We also show the zeros for the strips S^^m (black □) and Sz^q (red o). We depict the 
isolated limiting points with the symbol x. The zeroth-, first-, second- and third-order 
asymptotic expansions for the outward branches [cf. (7.14)] are depicted as dashed orange, 
dotted red, dot-dashed magenta, and dot-dot-dot-dashed blue curves, respectively. Please 
note that for m = 3 and k = 2,5, the first- and second-order asymptotics coincide, as the 
term of order has the factor sin[(2/i: — l)7r/3] = 0, which vanishes for k = 2,5. 
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Zeros sq lattice m = 4^,+^ ^eros sq lattice m = 4, 




Figure 12: Limiting curves for square-lattice strips of width m = 4 with two extra sites. 
We also show the zeros for the strips S'4,40 (black □) and (red o). We depict the 
isolated limiting points with the symbol x. The zeroth-, first-, second- and third-order 
asymptotic expansions for the outward branches [cf. (7.14)] are depicted as dashed orange, 
dotted red, dot-dashed magenta, and dot-dot-dot-dashed blue curves, respectively. 
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Zeros sq lattice m = 5p^^ 
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Zeros sq lattice m = 5 
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Figure 13: Limiting curves for square-lattice strips of width m = 5 with two extra sites. 
We also show the zeros for the strips S'5,50 (black □) and S'5,100 (red o). We depict the 
isolated limiting points with the symbol x. The zeroth-, first-, second- and third-order 
asymptotic expansions for the outward branches [cf. (7.14)] are depicted as dashed orange, 
dotted red, dot-dashed magenta, and dot-dot-dot-dashed blue curves, respectively. Please 
note that for m = 5 the zeroth- and first-order asymptotics coincide, as the term of order 
has the factor /i(5) = 0. 
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Figure 14: Limiting curves for square-lattice strips of width m = 6 with two extra sites. 
We also show the zeros for the strips Sq^q (black □) and 56,120 (red o). We depict the 
isolated limiting points with the symbol x. The zeroth-, first-, second- and third-order 
asymptotic expansions for the outward branches [cf. (7.14)] are depicted as dashed orange, 
dotted red, dot-dashed magenta, and dot-dot-dot-dashed blue curves, respectively. 



110 



Zeros sq lattice F++ bc's 2eros sq lattice F++ bc's 




Figure 15: Limiting curves Bm for square-lattice strips Sm,n of widths m = 1 (black), 
m = 2 (red), m = 3 (green), m = 4 (blue), m = 5 (pink) and m = 6 (orange) with two 
extra sites. The solid dark gray curve corresponds to the limiting curve for a strip of 
width m = 11 with cylindrical boundary conditions [34]. 
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Figure 16: Real part of the free energy Re fmiq) [cf. (5.2)] for square-lattice strips Sm,n 
of width m (and length n — )■ oo) as a function of Reg, for (a) Img = 3 and (b) Img = 4. 
Widths are m = 1 (black), m = 2 (red), m = 3 (green), m = 4 (blue), m = 5 (pink) and 
m = 6 (orange). The solid dots show the points of discontinuity in the derivative of the 
free energy. 
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Sq lattice F++ Im(q) = 0.5 

CO I 1 1 1 1 1 1 1 1 1 1 1 1 r 




Sq lattice F++ lm(q)=l 
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Figure 17: Real part of the free energy Refm{q) [cf. (5.2)] for square-lattice strips Sm,n 
of width m (and length n — )■ oo) as a function of Reg, for (a) Img = 0, (b) Imq = 0.25, 
(c) Img = 0.5 and (d) Img = 1. Widths are m = 1 (black), m = 2 (red), m = 3 (green), 
m = 4 (blue), m = 5 (pink) and m = 6 (orange). The solid dots show the points of 
discontinuity in the derivative of the free energy. 
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(a) (b) 

Figure 18: Vertex labels used in the proof of Proposition B.2, for (a) m even and (b) 
m odd. Vertices with a prime are obtained from the corresponding unprimed vertices 
by reflection Tl2. When m = 2£ + 1 is odd, the extra vertex (which is invariant under 
reflection) is denoted i + 1. 



114 



